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Abstract 



We present an effective Lagrangian formalism for the calculation of flavour 
changing neutral and charged scalar currents in weak decays including SU (2) x U (1) 
■ symmetry breaking effects and the effects of the electroweak couplings g\ and g2- 

We apply this formalism to the MSSM with large tan/? with the CKM matrix as 
the only source of flavour violation, heavy supersymmetric particles and light Higgs 
bosons. We give analytic formulae for the neutral and charged Higgs boson cou- 
plings to quarks including large tan/3 resummed corrections in the SU(2) x U(l) 
limit and demonstrate that these formulae can only be used for a semi-quantitative 
analysis. In particular they overestimate the effects of large tan (5 resummed correc- 
tions. We give also improved analytic formulae that reproduce the numerical results 
of the full approach within 5— 10%. We present for the first time the predictions for 
the branching ratios B® d ^> /U + /i~ and the B ds — B d s mass differences AM^ that 
include simultaneously the resummed large tan/3 corrections, SU(2) x U(l) break- 
ing effects and the effects of the electroweak couplings. We perform an anatomy of 
the correlation between the increase of the rates of the decays B® d — ► n + fi~~ and the 
suppression of AM S , that for large tan (3 are caused by the enhanced flavour chang- 
ing neutral Higgs couplings to down quarks. We take into account the constraint 
from B — > X s j clarifying some points in the calculation of the large tan f3 enhanced 
corrections to this decay. 



1 Introduction 



Supersymmetric theories are among the leading and most thoroughly investigated candi- 
dates for extensions of the Standard Model (SM). Even if supersymmetric particles are so 
heavy that they can be produced only in very high energy collisions to be studied at Teva- 
tron, LHC and NLC, their virtual effects can be investigated also in low energy processes 
like particle-antiparticle mixing, CP violation and rare decays of hadrons and leptons. A 
number of laboratories in Europe, USA and Japan will contribute to this enterprise in an 
important manner in this decade. 

Over the last fifteen years a vast number of papers discussing supersymmetric effects 
in such low energy processes has been published [1, 2, 3]. While K°-K° mixing, B®-B® 
mixing, CP violation in K — > tttt decays and in particular the radiative decay B — > X s 7 
played the leading role in these analyses so far, in view of forthcoming experiments an 
even more important role will be played in the future by the rare decays B — > X s e + e~ , 
B Q sd ^ n + n~ and K — > irvv, the B Q S -B Q S mixing and various CP violating transitions and 
asymmetries. 

In the version of the supersymmetric theory, in which the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix remains the only source of flavour and CP violation in the quark sector 
and tan (3 (the ratio of the two vacuum expectation values v u /vd) is between 1 2 and 10, the 
presence of supersymmetry does not drastically change the SM predictions for processes in 
question. While effects even as high as 30% could still be possible in AF = 2 transitions, 
generally smaller effects are predicted in rare K- and B-meson decays. Moreover, for su- 
persymmetric particles heavier than ~ 500 GeV their effects become unmeasurably small. 
Most recent analysis of this scenario has been presented in [4] , where further references to 
a rich literature can be found. We will call this scenario minimal flavour violation (MFV) 
at low tan /3, deferring the precise definition of our scenario to section 3. 

Much larger effects are still possible in supersymmetric models with new flavour vi- 
olating interactions originating from the misalignment of quark and squark mass matri- 
ces [5, 1, 6, 7, 8]. However, such models contain many free parameters and their predictive 
power is small. The same is also true for models with broken R-parity. 

In this paper we investigate another very interesting scenario, namely MFV in SUSY 
with large tan (5. It is theoretically appealing (e.g. it is consistent with the approximate 
unification of top and bottom Yukawa couplings at high energies predicted by some 5*0(10) 
models) and as predictive as MFV SUSY with low tan (3 (the number of free parameters is 
the same in both cases). In this scenario very large deviations from the SM predictions in 

1 For top squarks below 1 TeV tan/3 ^ 2 is most likely excluded by the unsuccessful Higgs boson search 
at LEP. 
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certain low energy processes are possible even for heavy supersymmetric particles but the 
calculations are more technically involved due to the need for resummation of large tan f3 
enhanced contributions. One of the aims of the present work is to formulate a general 
framework allowing to handle these complications. 

First discussions of large tan (3 scenarios in multi-Higgs models and supergravity mod- 
els can be found in [9] and [10], respectively. The importance of supersymmetric large 
tan (3 effects in processes involving charged Higgs particles H ± and flavour conserving 
processes involving neutral scalars has been also known for some time [11, 12, 13, 14, 15]. 
They originate in the tan (3 enhanced modifications of the standard relations between the 
original Lagrangian down-quark mass parameters m d , m s , m h (which determine the rel- 
evant Yukawa couplings) and the running ("measured") quark masses ma, m s , nib. For 
example, in the case of the 6-quark one has 

™b . , 

1 + e b tan (3 

The correction e& results from supersymmetric QCD and electroweak one-loop corrections. 
As |e&| can be of order 0.01, substantial enhancements or suppressions (depending on the 
sign of eb) of the Yukawa coupling are possible for tan/5 = (9(50). It has also been 
demonstrated [17] that expressing mj through nib by means of eq. (1.1) in the tree level 
couplings iH + b, t\ + b and in the flavour conserving neutral Higgs couplings resumms for 
large values of tan (3 dominant supersymmetric corrections to all orders of perturbation 
theory. Such a resummation is necessary for obtaining reliable predictions for measurable 
quantities [14, 16]. 

Even more important for phenomenology is the fact that for large values of tan/3 
flavour changing neutral currents in the down-quark sector mediated by Higgs scalars 
can be significantly enhanced. In the SM such currents are also induced at one loop by 
tW ± penguin-like diagrams with Z° replaced by the neutral Higgs boson but contribute 
negligibly to the rare processes. In supersymmetry with large tan (3 neutral Higgs boson 
penguin-like diagrams with charginos and stop-quarks in the loop can be more important 
than the standard Z Q penguins thanks to tan (3 enhancement of the down-quark Yukawa 
couplings. They have been first considered in [18] in connection with the B — > X s l + l~ 
decay and in [19] and subsequently found to increase by orders of magnitude the branching 
ratios of the rare decays B° s d — > [20, 6, 21, 22] and to decrease significantly the B° s -B° s 

mass difference AM S [23] relative to the expectations based on the SM. 

Two additional related aspects of large supersymmetric corrections in the down-quark 
sector are the following. Elements of the physical CKM matrix, to be called Vjf in what 
follows, differ from Vji present in the original Lagrangian by calculable flavour dependent 
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corrections fji [24, 19, 20]: 



Vji = Vjffjj (1.2) 

(J, / denote flavour indices), fjj can substantially deviate from unity for large tan/3 and 
can in principle affect all processes in which the CKM matrix is relevant. 

Finally, as has been shown in [25] in the context of the analysis of B — > X s 7 decay, 
there are additional large tan /3 enhanced corrections in the charged Higgs couplings, that 
in ref. [25] have been parametrized by e' b and e' t . 

It is desirable to calculate all the four effects, that is: 

1) Large tan/3 effects related to (1.1), 

2) Enhanced neutral scalar-penguin diagrams, 

3) Large tan/3 effects related to (1.2), 

4) Enhanced corrections to charged Higgs vertex diagrams 

in a self-consistent framework and to include them in a phenomenological analysis cor- 
relating predictions for B — > X s j, B° sd — > rates and for B° sd -B° sd mass differences 
AM s>d . It is also important to formulate general rules allowing to include large tan/3 
effects that can be used in other processes not considered in our paper. 

Several steps towards this goal have been already made. In refs. [17, 25, 26] large 
tan/3 effects related to 1) and 4) have been discussed in the context of the B — > X s ^ 
decay and recipes for including them in the charged Higgs (H^) and charged Goldstone 
boson (G^) vertices have been formulated for specific quark flavours. In another elegant 
analysis [27] the effects l)-3) in the case of B° s d — > fi + decays and B° s d -B° s d mixing 
have been calculated in the SU(2) x U(l) symmetry limit, neglecting the electroweak 
couplings #2 and g\ in comparison with a s , y t and y b . These authors confirmed sizeable 
enhancements of BR(B° s d — > yU + /i~) and the suppression of B° s -B° s mixing pointed out 
previously in refs. [20, 6, 21, 22] and [23], respectively. 

During the final stages of completion of our paper a model independent analysis of 
rare processes in theories with the CKM matrix as the unique source of flavour and CP 
violation has been presented in [28]. While those authors also investigated large tan/3 
effects in BR(B° sd — > fi + fi~), AM S and B — > X s 7, their analysis was performed in the 
same SU{2) x U{\) symmetry limit assuming the dominance of a s , y t and 

In the present paper we go beyond these analyses by 

• Calculating all the four effects in an effective Lagrangian approach that goes beyond 
the SU(2) x U(l) symmetry limit considered in [27, 28] and includes the effects of 
the electroweak gauge couplings. 
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• Calculating all the four effects in the SU(2) x C/(l) symmetry limit with vanishing 
electroweak gauge couplings thereby confirming and, in certain cases, correcting and 
generalizing the analytical rules for the inclusion of large tan (3 effects presented in 
refs. [25, 26]. Our results agree with those of [27, 28] if the factors denoted here by 
€q, ey, €q and e' Y are assumed to be flavour independent. 

• Analysing numerically the validity of the above approximation and giving simple 
recipes which improve it significantly. 

• Including all these effects in an analysis of the branching ratios BR(B ds — > 
B®-Bg mass difference AM S and of correlations between them taking into account 
the constraint from B — > X s 7. 

To our knowledge no complete analysis of all these topics has been presented in the 
literature so far. While the recent analysis in [28] considered many of these issues and 
went beyond supersymmetry, it was performed in the 577(2) x U(l) symmetry limit with 
9i — 92 = 0, the approximation the accuracy of which we want to investigate here. In 
particular it is interesting to see how big effects can still be expected in B® s fi + fi~ 
given the constraints on the magnitude of the scalar flavour changing neutral currents 
stemming from the B®-B® mixing [23]. This should also allow to assess the claim made 
on the basis of approximate formulae by the authors of ref. [27] that the effects found 
previously in [6, 22] and [23] can be decreased by a factor as large as 5 when the large 
tan/3 effects described by (1.1) and (1.2) are properly taken into account. 

In a recent letter [29] we have presented some of the formulae resulting from our 
analysis and in particular we have analyzed the correlation between the increase of the 
rates of the decays B® d — > ji + fx~ and the suppression of AM S . In this paper we present 
the details of our formalism, derive all results and extend the phenomenological analysis. 

Our paper is organized as follows. In section 2 we present general formalism for 
calculating the Higgs flavour changing and flavour diagonal vertices based on the effective 
Lagrangian approach valid beyond the SU (2) x U (1) symmetry limit. In section 3 we relate 
our approach to the one of [27, 28] and present explicit expressions for flavour violating 
neutral scalar (S°) and charged scalar (H ± , G^) couplings in the SU(2) x U(l) symmetry 
limit. Subsequently we generalize these expressions to allow for flavour dependence in the 
parameters e , ey, e' and e' Y . We compare our formulae with those present in the literature 
in section 4. In section 5 we asses numerically the importance of the corrections calculated 
in sections 2 and 3. We also investigate the effects of SU(2) x U{\) breaking and of the 
effects of non-vanishing electroweak couplings comparing our results with those present 
in the literature. This includes in particular the parameters e and e' of refs. [25, 26]. 
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Finally using the flavour dependent formulation of section 3 we find analytic formulae 
for various couplings that within 5 — 10% reproduce the numerical results of the full 
approach presented in section 2. In section 6 we collect the basic formulae for AM^ S , 
BR(B® S — > and BR(B — > X s ^) in supersymmetry at large tan/3. Section 7 is 

devoted to the numerical analysis of these quantities and the investigation of possible 
correlations mentioned above. We conclude in section 8. 

2 General formalism 
2.1 Effective Lagrangian 

The most efficient way of handling heavy particle effects in low energy processes involving 
quarks, is the effective Lagrangian technique. Let us assume that a sector of heavy 
fields in a theory (e.g. sfermions, gluinos, charginos and neutralinos in the MSSM) can 
be integrated out without violating gauge invariance. We consider processes occurring 
at energies low enough so that the amplitudes cannot have imaginary parts related by 
unitarity to particles belonging to the heavy sector. Renormalized in the MS scheme 
with some renormalization scale Q, virtual effects of the decoupled particles (Q should be 
chosen to be of the order of their masses) corresponding to self-energy diagrams in Fig. 1 
can be summarized by the effective kinetic and mass terms of light fermions: 

- (dj) R (m d + £™l(0)) (d/) L - {dj) L (m d + E^(0)) (d T ) R (2.1) 

plus similar terms for the up-type quarks and leptons. Here I, J are flavour indices with 
di = d, d 2 = s, d 3 = b. Analogous notation will be used for the up-type quarks. Exact 
expressions for quark self-energies are listed in the Appendix A. 2. 

In order to make the formulae more compact we denote (in some of them) the diagonal 
quark mass matrices simply by = diag(m ( ; 1 , m<2 2 , "Ms) = diag(mrf, m s , mj,), m u = 
diag(m ui , m U2 , m U3 ) = diag(m M , m c , m t ). Hermiticity of the effective Lagrangian ensures 
that (E^)t = j« Li (E ^)t = j« R and (S ^)t = 

Similarly, virtual effects of the decoupled particles corresponding to fermion-gauge 
boson vertices of the type shown schematically in Fig. 2 can be summarised as 

4ff = -IWlI" (Ft" + ^F[ d ) JI {d:) L Zl - Jdj)- Rl » (Fi d + AF z R d ) U (dj) R Z° 

- Mil 1 * (Fl u + AFf «) JI ( Ul ) L Zl - J^f (Fi u + AFi u ) 31 { Ul ) R Zl (2.2) 

- MlY {FY + AFf ) JI {d:) L W; - J^jfrf (Fg + AF™)" ( dj ) flW J + H.c. 
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-i (E q VL i>P L + E% i>P R 



Figure 1: One- loop threshold corrections to fermion propagators. Pl = Pr = 
and q = d or u. 

where (F^ q ) JI oc 5 JI , (F^) J/ = {g/y/2)Vn and, of course, Fjf = 0. 




-i>f (AF/ 7 P L + AF^Pfl) 



Figure 2: One-loop threshold corrections to fermion-gauge boson vertices. 

Finally, other virtual effects of the decoupled heavy particles can be represented in the 
effective Lagrangian by various effective vertices which are operators of dimension higher 
than four (non-renormalizable operators) such as four-fermion operators or dimension 5 
(chromo)magnetic operators etc. 

2.2 Diagonalization 

The next step is to rescale the quark fields so to render their kinetic terms canonical: 

+ \dj) L , ( dl ) R ^(l + h: d v ^ Ij \dj) R , (2.3) 

and the same for the up-type quarks. After this operation the quark mass terms become 
£ mass = -(dj) R (m d + Am d ) JI (d^L - (dj) L (m d + AmJ) (dj) R (2.4) 

where 

Am, ee S^ L (0) + ^ d VR (0)m d + ~m d S^(0) . (2.5) 



Similar expressions are obtained for the up-type quarks. 

The mass terms (2.4) can be diagonalized by two independent rotations: 



d L -> B L d L , d R -> B R d R (2.6) 

so that 

'T>l(m d + Am d )D L ] = m dj 5 JI (2.7) 

with m dj denoting the corrected mass eigenvalues. To one-loop accuracy the transforma- 
tions (2.6) read: 

d L ^(l + AD L )4 d R ^(l + AB R )d R . (2.8) 

The unitarity of ^l(r) requires that (up to higher order terms) ADL™ = -AD L ( fi j and 
AD^ fi j = 0. In this approximation the corrected mass eigenvalues are given by 2 

m dj = m dj (1 + K dj ) . (2.9) 

where 

Kdj = ~^ + 9 VR + 9 VL ^ = ~^ • ^ °> 

lll dj Z, Z, lll dj 

Formula (2.9) generalizes (1.1) to arbitrary down-quark flavours. 

In order to find the matrices AD L( - R ) we follow ref. [27] and decompose Am d explicitly 
into its diagonal and off-diagonal parts 

(Am/ 7 = m djKdj 5 JI + (A'm d ) JI (2.11) 

where by definition (A'm d ) JI = (Am d ) JI for J ^ I and (Am d ) JJ = 0. 
The condition determining the AD^^) takes now the form 

(AD^ + IrfADi + Am d ) JI = (2. 12) 

for J I. Note that following ref. [27] we multiply by AD^ also the tan/3 enhanced 
diagonal correction to the down-type quark masses, i.e. include terms which are formally 
of higher order. The non-zero entries of AD[j fi j are then 

ADf = .M< + (Amt)"m d , (213) 

m dj m di 

2 If the diagonal entries of the matrices are complex, the rotations (2.8) have to be supple- 

mented by additional chiral rotations of the quark fields in order to get real and positive quark masses. 
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m dj - m dj 



Using the hierarchy of quark masses: m da ^> m d2 ^> m dl (m,d 3 3> ^ m dJ and the 
fact that Amf oc (AroJ) J/ oc m dl (see section 3 for explicit expressions) we find for 



J > I 



ADf = - J-Arn;" + f ^) , 

AD- = -^(A^)" - |Amf ~ O (2.15) 



and for J < 7 

ADf = J-(AmS)'" + o( r ^V 

AD- = lAmf + 5(A™1) J/ - ^ (^) (2-16) 

Note that AD L is 0(1) with respect to the masses but AT) R is always suppressed. This 
suppression can be occasionally removed by a large mass multiplying AD R . 

The same formulae (with appropriate replacements d — > u, D — > U) hold also for the 
up-type quarks. In particular we have 

l4(m u + Am u )U L ] JI = m Uj 5 JI (2.17) 

with m Uj denoting the corrected mass eigenvalues. However the correction Am u is not 
enhanced for large tan (3 (see the discussion following (3.16)) in contrast to Am d that is 
enhanced in this limit. Consequently one can set m Uj ~ m Uj and the condition analogous 
to (2.12) reads 

(AU t R m u + m u AU L + A'm M ) J/ = (2.18) 

for J ^ I. The explicit expressions for AJJ^^ are then obtained directly from (2.13)- 
(2.16) by making appropriate replacements (jn dj — > m Uj , etc). 

Obviously, these are the masses m dj , and m Uj m m Uj that have to be identified 
with the running mass parameters of the low energy effective theory valid for scales 
~ Mz- Parameters m dj , fn Uj are therefore directly related through the known QCD 
renormalization group running to the low energy quark mass parameters. 
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2.3 Effective CKM matrix 



Rotations (2.8) have the effect of renormalizing the CKM matrix. After the operations 
(2.3) and (2.8) the W-boson vertices become 

= -M^ (*f + AFf ) J/ (d 7 ) L ^; - (AFjT) " (d,)*^ + H.c. (2.19) 

with 

AFf = AFf + \j% L F? + ±F?1$ L , AFg = AFJ (2.20) 

and 

(^f) J/ = fvft - ^ (XJtVD £ ) U « ^ (V + AV[V + yAD L ) J/ . (2.21) 

Since the corrections AF^ R are negligible (for example, for all JI the gluino exchange 
computed in [33] gives AFf /Ff < 3 x 10~ 3 ), it is the matrix Vjf that has to be identified 
with the CKM matrix whose elements are determined from the low energy processes. The 
matrix Vji is the CKM matrix of the MSSM which for a given spectrum of sparticles is 
related in a calculable way to the measured matrix Vjf . Thus, the formula (2.21) is the 
realization of the relation (1.2). The correct procedure in phenomeno logical applications 
is then to calculate first all relevant quantities in terms of Vjj and subsequently express 
the latter in terms of Vjf that should be determined from the measured branching ratios. 

It is known [24, 20], that for large tan/5 the ratios of the elements V£f /V$i ~ Vff /Vis 
can, for 1^3, substantially deviate from unity. Simple approximate (but very accurate) 
expressions for fjj = Vji/Vjf defined in eq. (1.2) will be given in section 3 (see eq. 
(3.53)) and the corresponding numerical analysis in section 5. 

The elements |V^&| and \V c b\, that are affected by these corrections are usually deter- 
mined from tree level decays under the assumption that new physics contributions to the 
relevant branching ratios can be neglected. This assumption is clearly violated in the 
case of supersymmetry at large tan/3. However, in this case the most important tan/3 
enhanced loop corrections to tree level decays, can be absorbed in the V cS . Therefore, 
the quantities which experimentalists extract from these processes can be identified with 
|Kf | and \Vf\ defined in (2.21). 

On the other hand, as we will see the ratio is to an excellent approxima- 

tion not affected by these corrections and |Kib|/|Vd>| ~ iKf^l/lK^I- ^ should also be 
emphasized that 

V cS (V cS Y = 1 (2.22) 

so that all unitarity relations, in particular the formulae related to the unitarity triangle, 
expressed usually in terms of Vjj, are also valid for Vjf. 
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For the Z° couplings the rescalings (2.3) amount to the replacements 

lAr L(R) ~^ ^ r L(R) — L(R) 2 VL ( R ) L ( R ) 2 VL(R)- l Z ' ZO J 

To one-loop accuracy, the subsequent unitary rotations (2.8) do not affect the Z° vertices. 
Flavour changing introduced in higher orders can be treated along the lines discussed in 
subsection 2.6. 

2.4 Neutral Higgs-fermion vertices 

The one-loop corrections to scalar (S°) neutral Higgs boson vertices induced by the de- 
coupled particles can be written as 

£5 = -¥jVr (F? + AFf ) JI ( dl ) L S° - (ff + AFf ) J/ (rf 7 )^ (2.24) 

and similarly for the up-type quarks (see Appendix A. 3 for the expressions for AFl R ). 
For neutral Higgs and Goldstone bosons we have 




Figure 3: One- loop corrections to neutral Higgs-quark vertices. 



(F d L f R) ) JI = 4\ R) m dj 5 JI ee (F d L f R) Y5 JI . (2.25) 

with explicit expressions for rfjj, R \ given in (3.20) and (3.21). After the operations (2.3) 
and (2.8) the expressions inside the brackets in eq. (2.24) become respectively 

pdS + AF dS + l^d^dScorr + I^corryjd ^ + ^t^corr + ^Scor^j^ 

F R + AFf + ^ L Ff corr + ^Ff corr £^ R + AD[Ff corr + F R Scovv AD R (2.27) 
where we have defined 

K«) rr ) J - 4 S iR) m dj + (AF£f R) ) JJ . (2.28) 
10 



Of course, in the phenomenological applications one has to express Fjf in terms of the 
"observable" by using eq. (2.9): 



j 

ji 



(2.29) 



with an analogous formula for (Fr S ^J . In writing the second expression in (2.29) we 
have assumed that the vector parts of the fermion self energies do not contain large 
contributions and neglected terms of order e d (pvL + ^vr) an d smaller. We have also 
defined the quantity e dj 



. \ J J 



e dj = V " %L ' R (2.30) 
m dj tan p 

which for J = 3 corresponds to ej, of ref. [25]. It is also instructive to make contact 
with the standard perturbative (on-shell) calculation of the one loop dS°d vertex which 
gives [6, 22] 

V fm dj S JI + (AFt s ) JI -v d L S ^l. (2.31) 

By expanding (2.26) strictly to one loop (neglecting the difference between m d and Tn d 
in all 1-loop terms) and using the explicit expressions (2.13), (2.14) as well as (2.9) and 
(2.10) it is easy to see that the two approaches indeed coincide. In particular, to 1-loop 
accuracy the vertex (2.26) is independent of the vector self energies, both for the diagonal 
and off-diagonal transitions. 

On the other hand, keeping e dj in the denominator in (2.29) resumms dominant cor- 
rections of order a s (/i/mg) tan/3 to all orders of the perturbation expansion [17] which are 
important for neutral Higgs decays into bb or t + t~ pairs [13, 14, 16]. Thus, the diagonal 
couplings of scalars to the down-type quarks are 




Ctr = -(dj) R rjf - ■ d -[ - + AFf (dj) L S -R.c. (2.32) 



The expression for the flavour non-diagonal vertex can be simplified by taking into 
account the hierarchy of masses m dj : m dz ^> m d , 2 ^> m dl (rn d3 ^> m d2 ^> m dl ) and using 
the expressions (2.15) and (2.16) for the AD L ( R ) matrices. One has also to take into 
account that (as we will see shortly) in the limit of SU (2) x U(l) symmetry (Am d ) JI oc m dj 
and (Am d ) JI oc m dl . In general one gets 



£ off-dia g = _ {dj)R r x 5l JI {di)L gO _ {dj)L 1 X S 1 - {di)RS (2>33) 



JI 
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where for J > I 



X RL 



.11 



( AFf + -Zi R F? Scon + -F 



dScorrs^d 



.11 



^pdScoiT^ 1 

m dj 



(Am d ) JI (2.34) 



X 



LR 



(^jpdS I \_Y^ pdScorr . pdScorr j-i 
i? 2 ^ 2 ^ 



(2.35) 



ft 



c(Scorr\ ^ 



/ ZT'dScorrA ^ / ZT'dS'corrA 



(Ami) 



and for J < I the expressions for (X RL ) JI and (Xf R ) J/ can be obtained by using the rules 
(remember that (Ff corr ) J * = (Ff corr ) J ): 



X 



RL 



JI 



X 



LR 



.11 



X 



LR 



IJ* 



X 



LR 



JI 



X 



St 
RL 



.11 



X 



RL 



IJ* 



(2.36) 



Explicitly for J < I one gets 



X RL 



JI 



(af? s 



1. 



1 

ft 



dScorr _|_ ^ pdScorr ^d 



VL 



.11 



(pdScorA 1 ( f F dScovv\ J ( F dScorry m 



(2.37) 



UJi 



m dl 



m dl 



di 



X 



LR 



1 1 

— d rpdScorr , rpdScorrs^d 

2~VL P R "+" 2 R yi? 



-^ r L(AmS)". (2-38) 



We observe that for J > I (I > J) the effects of the rediagonalization in Xlr {Xrl) 
are suppressed by a factor ~ m dl j 'm dj (m dj j 'm dl ) compared to its effects in X RL (X LR ). 
It has to be remembered that m d and V entering the calculation of AFf^ R y Am d etc. 
have to be ultimately expressed in terms of m d and V cS . 



2.5 Charged Higgs-fermion vertices 

In the same manner the operations (2.3) and (2.8) imply the following charged scalar 
(H^) and charged Goldstone boson (G^) couplings 
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Figure 4: One-loop corrections to charged Higgs-quark vertices. 



where (i = H+, G+) 
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Here 



AF k L 
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r] R Vm d 



(2.41) 
(2.42) 

(2.43) 
(2.44) 

(2.45) 



and AF*™ defined in analogy with (2.2) (see Appendix A. 3 for the explicit expressions). 
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(2.46) 



(2.47) 



In writing (2.41) and (2.42) we have used (2.21) in order to convert the CKM matrix 
in the tree level vertices Fhm the effective CKM matrix V eS . 

In what follows it will be useful to define the parameters ejj , e^j R , Ejf and Cjf" through 



pH 

RL 



P 



G 
RL 



.11 



.11 



vZm UJ V]f(i-e^ 



GL\ 



pH 
LR 



LR 



.11 



.11 



riZVJfm^l-e™) 



GR\ 



(2.48) 
(2.49) 



This parametrization of the H + and G + couplings differs from the parametrization 
in terms of e' factors introduced in ref. [25] but is more natural for our purposes. In 
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particular it emphasizes the fact that ultimately all vertices should be expressed in terms 
of rUi and Vjf. 

The numerical analysis of the parameters €jf, €jP, e'jf and e^P is presented in sec- 
tion 5. While the general formulae for these parameters are rather involved, they simplify 
considerably in the SU{2) x U{1) symmetry limit with the gauge couplings g\ and #2 set 
to zero as we will see in section 3. In particular the parameters ej/ and ej/ 2 vanish in 
this limit. Therefore ej/ and e^f should be generally smaller than ejf, and €jj R . This 
expectation is confirmed by our numerical analysis presented in section 5. 



2.6 Chargino-Quark-Squark Vertices 

In the approach formulated above charginos and squarks are integrated out and do not 
appear as explicit degrees of freedom in the low energy effective theory. Their virtual 
effects are accounted for by corrections Am djU , various formfactors AF and, as already 
mentioned, by higher dimension operators. Let us consider for definiteness one of several 
such non-renormalizable operators 

A£ cff = _ C VLR n (jdjj2^( dl ) L ) (jdMha^R) (2.50) 

which is generated by chargino-stop box diagrams in the course of integrating out these 
particles. Cjf^ are the Wilson coefficients that depend on flavour indices J, J, M, N and 
can be calculated by using the original MSSM vertices 

£ mssm = _ £ — ~ t + + h.c. (2.51) 

i=l,2 

which (in the matrix notation) have the following structure [30]: 

4= + ~Vv, c R = — .V.m d (2.52) 
\ sin p J cos p 

Here a 1 , b l and h % are 6x3 matrices that depend on chargino compositions and left-right 
mixing of squarks. The coefficients like C V jl^ N corresponding to definite chiralities of 
external quark fields are given as sums over squark and chargino indices of appropriate 
products of d L , c R , (c l L y and (c R Y weighted by some functions of chargino and up-type 
squark masses. 

The operations (2.3) and (2.8) performed on the quark fields in the effective Lagrangian 
affect also the coefficient C^f^N- Neglecting rescalings (2.3) which are not enhanced by 
tan/3 it is easy to see that the rotations (2.8) amount simply to the multiplication from 
the right of each c l L and c l R entering CJj L m N by the matrix D L and D R) respectively. From 
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(2.21) it follows that in the approximation AUl(r) ~ this has the effect of replacing V 
by V cS in all c l L . For the factors c l R the effect of the rotation is more complicated: 

V ■ m d ■ D R « V cS ■ D{ ■ m d ■ D R « \/ cff • m d + \/ cff • (m d • AD, - AD L ■ m d ) 

Thus, the most important effects of the rotations (2.8) in the non-renormalizable 
operators can be simply taken into account by calculating their Wilson coefficients using 
the modified chargino-squark couplings: 

^{* + » (2 ' 53) 



V cS -m d + V cS -(m d - AD R - AD L -m d )) (2.54) 



^ R cos/3 

where m Uj w m Uj and m dl = m dl /{l + e dl tan 

The effective vertices of the effective Lagrangian calculated as described in this section 
do depend on the renormalization scale /j,r. According to the standard reasoning this 
scale should be of the order of the mass scale of the particles which are integrated out. 
The couplings and masses of the effective Lagrangian should be then treated as running 
parameters. For example, the contribution of the Higgs boson exchanges to the Wilson 
coefficients of the effective Hamiltonians discussed in sec. 6 should be calculated at the 
scale uri ~ Mffiggs ~ ■ If the MgusY ^ ^Higgs ~ Mh+ , the couplings of the effective 
Lagrangian should be evolved with the effective theory RGEs from /j, R ~ Msusy down 
Uri ~ Mniggs ~ Mh+. In our numerical results we neglect this running as it would not 
change the results drastically. It should be however taken into account to predict physical 
quantities in the MSSM with higher accuracy. 

Finally let us stress that although in what follows we will restrict our attention to the 
version of the MSSM in which flavour and CP violation are ruled by the CKM matrix, 
the effective Lagrangian approach described in this section is more general and can be 
used [6, 7] also in the case of flavour and CP violation originating in the squark sector. 



3 The limit of unbroken SU(2) X U(l) symmetry 

In this section we will discuss the approach of ref. [27] that can be considered as a special 
limit of the more general approach presented in section 2. We will also compare our 
results with the recent analysis presented in [28] that goes beyond supersymmetry but 
similarly to [27] is based on the same approximation. In order to relate these approaches 
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to our approach we consider decoupling of sparticles in the limit of unbroken SU (2) x U(l) 
symmetry. The electroweak symmetry breaking is then taken into account after sparticles 
are integrated out. The main goal of this section is the derivation of explicit expressions 



for the flavour changing neutral scalar couplings 



scalar couplings 
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r RL 



Jl 
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r LR 



Jl 



pG 
r RL 



Jl 



X RL 



and 



pG 
r LR 



Jl 



Jl 



and 



X LR 



ji 



and the charged 



at large tan (3 in the SU{2) x U(l) 



symmetry limit. This approximation should be valid if the sparticle mass scale is larger 
than the mass scale of the Higgs boson sector (set by Mh+). It should be emphasized 
that the absence of vacuum expectation values before decoupling implies neglecting the 
left-right mixing of squarks even for non- vanishing A u ^ and/or \x parameters. 

In this section following the practice in the existing literature we will assume also the 
dominance of a s and top and bottom Yukawa couplings neglecting the contributions of 
the electroweak couplings g-i and g±. We will assess the validity of this approximation in 
section 5. 



3.1 Minimal Flavour Violation and Flavour Dependence 

The large tan (3 enhanced corrections in the approach of [27, 28] are governed in the MSSM 
by four universal parameters that we will denote by eo, ey, e' and e' Y - These parameters 
as defined in [27, 28] are flavour independent. Therefore flavour dependence enters the 
Higgs and chargino couplings only through the CKM matrix V e g. This scenario can only 
be realized if the soft SUSY breaking mass squared matrix 171% is strictly proportional 
to the unit matrix. This corresponds to the Minimal Flavour Violation as defined in ref. 
[28] and to the scenario A in [22]. 

However, in more realistic situations the diagonal entries «^q 3 , m lj 3 and m 2 D3 of the 
matrices rug, mfj and m 2 D that are related to third generation squarks differ from the 
entries corresponding to the first two generations to be denoted collectively by m 2 . As 
pointed out in [1] in this situation some flavour violation (still ruled by the CKM matrix) 
unavoidably appears in the up- or down-type (or in both) squark mass squared matrices. 
For m 2 7^ m Q 3 , at the phenomenological level, it becomes a matter of choice in which 
basis uiq is assumed to be diagonal. The two simplest choices are 

(M 2 D ) LL = m 2 Q + . . . and {M 2 v ) hh = V m 2 Q V^ + . . . (3.1) 

or 

(M 2 D ) LL = V^m 2 Q V + ... and (A4)ll = m 2 Q + . . . (3.2) 

where rriq is diagonal and the ellipses stand for other diagonal terms. In the classification 
in [22], the settings in (3.1) and (3.2) correspond to scenarios B and C, respectively. 
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In our calculations in the approach of section 2 we choose the soft SUSY breaking 
mass parameter m 2 Q as in (3.1). The scenario (3.2) with flavour violation in the down-type 
squark mass matrix would require the inclusion of box and Higgs penguin diagrams with 
gluinos and is beyond the scope of this paper. Note that in SUSY breaking scenarios like 
minimal SUGRA or gauge mediation, in which proportionality of itlq to the unit matrix 
is assumed at some high scale, the renormalization group evolution produces calculable 
off diagonal terms both in (A^^)ll and in (M 2 d )ll- This is however a secondary effect 
as the evolution has to generate first the splitting of the diagonal entries of the squark 
mass matrices. The magnitude of the off diagonal entries is then proportional to this 
splitting. Due to this fact the effects of the splitting are always more important in flavour 
changing amplitudes. Our framework should therefore constitute a good approximation 
to predictions of the SUSY breaking scenarios in which at some high scale the matrices 
uiq, mfj and m 2 D are proportional to the unit matrix and the trilinear soft parameters A t 
and Ab are real. 

Allowing for different values of the parameters mq 3 , m^, m 2 D3 and m 2 introduces some 
flavour dependence into eo, ey, e' , e' Y and the formulae for various couplings presented in 
[27, 28] have to be generalized to account for it. However, for the sake of comparison with 
[27, 28] we will for the most part of this section treat e , ey, e and e' Y as universal flavour 
independent quantities and only at the end, in subsection 3.6, we will provide formulae 
for various couplings that are also valid for flavour dependent eo, ey, e' and e' Y . As we 
will see in sec. 5, obtaining a numerically reliable approximation to the results of the 
complete calculation done as in sec. 2 will require extending the approach of refs. [27, 28] 
by taking into account both the flavour dependence of eo, ey, e and e' Y and the effects 
of the electroweak couplings gi and g\. The latter will be automatically included in the 
formulae presented in sec. 3.6 by defining appropriately flavour dependent factors eo, ey, 
e' and e' Y as discussed in section 5. Thus, the analytic expressions of sec. 3.6 are the basis 
of the approximation that within a few percent reproduces the numerical results of the 
full approach of section 2 which includes automatically both the SU(2) x U(l) breaking 
effects and the effects of electroweak couplings. 

3.2 Effective Lagrangians 

In this approach below the sparticle mass scale we get the effective Lagrangian whose 
part relevant for the neutral Higgs boson couplings to the down-type quarks has the 
form [31, 20] 

C cS = -eijH^dR ■ (Y d + A d Y d ) ■ q jL - H^*^ ■ A u Y d ■ q iL + H.c. (3.3) 
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where £21 = — 612 = 1 and are Yukawa coupling matrices. As in the papers [27, 28] 
we make further assumption (which by itself is not related to the symmetry limit) that 
the effects of the electroweak gauge couplings g± and #2 are negligible. We will examine 
numerically the validity of this approximation in sec. 5. With this assumption the only 
diagrams giving rise to the correction A u Y d are shown in Figs. 5a and 5b. The loop 
induced term — e^H^ A d Y d is always subleading in the large tan/3 limit and can be 
neglected. In the basis in which 

Y, = diag(y d ), Y u = diag(y u ) • V (3.4) 

the correction A u Y d can be easily seen to have the structure [27] 

(A u Y d ) JI = -y dj (e 5 JI + e Y y 2 t V 3J *V 31 ) = -y dj (ej5 JI + e y y t 2 A ") (3-5) 

where the contributions proportional to y\ and y 2 c have been neglected, 

lj = e + t Y y 2 t V 3J *V 3J » e + e Y y 2 5 J3 (3.6) 

and 

X Ji = V 3J *V 31 for J^I and A^ 7 = for J = I. (3.7) 

The approximation in (3.6) is justified by the hierarchical structure of the elements of 
the CKM matrix. The correction A d Y d has the same structure (3.5) as A u Y d but is 
subleading for large tan (3. 

As follows from figs. 5a and 5b, in (A u Y d ) JI the correction e depends on the gluino 
mass and the masses of the Dj and Qj squarks, while e Y depends on the \i parameter 
and the masses of the Qj and C/| squarks. Note that in this limit the parameter \i plays 
the role of the mass of the higgsinos. The dependence of e (Dj,Qj) and e Y (Q,j,U^) on 
the generation indices can be neglected if the squark masses are not significantly split. It 
will be taken into account in sec. 3.6 

Decomposing the two Higgs doublets as 

H[*> = -3 + ^ (c«H° - s a h° + i Sf3 A° - ic,G°) 

iff* = -3 + M° + cah° - icpA - lS ,G°) (3.8) 
(where c a = cos a, etc.) we find first 

A ( nL = d R ■ (^y d + ^j=A d Y d - ^uY d ) ■ d L + H.c. (3.9) 
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Figure 5: Vertex corrections in the SU(2) x U(l) symmetry limit for vanishing electroweak 
gauge couplings. Diagrams a) and b) give rise to corrections (A u Y d ) JI . Corrections 
(A ( iY c i) JI are generated by similar diagrams with outgoing replaced by incoming 
and factors /i* and A* u in diagrams a) and b) replaced by A d and fx, respectively. 
Diagrams c) and d) give rise to corrections (A d Y u ) JI . Corrections (A U Y U ) JI are generated 
by similar diagrams with outgoing replaced by incoming and factors fi* and A* d 
in diagrams c) and d) replaced by A u and /i, respectively. 
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Comparison with (2.4) gives then 

(Am d ) JI = (A d Y d - tan/? A u Y d ) J/ » -| tan/3 (A u Y d ) J1 . (3.10) 



In our conventions 

V2 Vdj ' muj ~ V2 
and v 2 d /cos 2 (3 = v 2 ^/ sin 2 (3 = l/y/2G F « (246 GeV) 2 . Consequently 



« U /O 1 1 \ 



(Am/ 7 = m dj tan/3 (ij5 JI + eyy 2 ^ 1 



(Ami) J/ = m d/ tan /3 (eV J/ + e Y |/ 2 A J/ ) (3. 12) 

K dj ~ e d,/ t an ft = lj tan /3 

with defined in (2.30) and e,j given in (3.6). The last equation tells us that in the 
approximation considered here e dj = ij. Thus using (2.9) we have 

™* = irfe- (3 ' 13) 

with m^j denoting the running quark masses as discussed in Section 2. Note that the 
large tan/3 corrections to m d and m s in contrast to m b do not involve the top Yukawa 
coupling. 

Similarly the corrections to the up-type quarks Yukawa coupling can be summarized 
by the effective Lagrangian 

C cS = -e l3 H\ u) u^i ■ (Y u + A U Y U ) ■ q jL - H\ d) *u^ ■ A d Y u ■ q lL + H.c. (3.14) 

The correction A^Y U is generated by the diagrams shown in figs. 5c and 5d. In the basis 
(3.4) one has 

(A,Y tt ) J/ = y Uj Vjj (e' + e' Y y 2 di ) . (3.15) 

where the factor e' Q depends on the masses of gluino and Uj and Qi squarks, while e' Y 
depends on the /i parameter and the masses of Qj and Dj squarks. For squark masses 
not significantly split the generation dependence of e' (Uj,Qi) and e' Y (Q i , Dj) on the 
generation indices can be dropped. The correction A U Y U has the same general structure. 
The mass term for the up-type quarks takes then the form 

C£L = -ur ■ {^VuV + ^uY u + ^j=A d Y^j ■ u L + H.c. (3.16) 

Note that the effects of the correction A d Y u are suppressed by tan (3 with respect to the 
effects of the correction A„ Y u which itself is not tan (3 enhanced with respect to the tree 
level mass term. The corrected mass term (3.16) is diagonalized by the rotations 

u L -> V ] ■ U L u L , u R -> u R U R (3.17) 



20 



that with Ul,r = 1 + AUl^ satisfy the conditions (2.17) and (2.18) with 
(Am u ) JI = ^= (A U Y U • + cot/? A d Y u ■ 



(3.18) 



Since Am u is not tan/5 enhanced, Ul,r are in general very close to unit matrices. Con- 
sequently 

m Uj m Uj . (3.19) 



3.3 Neutral Higgs-fermion vertices 

Using the effective Lagrangian (3.3) and the decomposition (3.8) it is easy to identify in 
the SU(2) x U(l) symmetry limit with zero g x and g 2 the vertex corrections AFffL for 
S° = H°, h°, A and G° defined in eq. (2.24). We find: 
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_ m dj T S 
Ji _ m dj 



Vd 



= {*R 

(ij5 JI + e Y y 2 X 



■2\-lf\ X S 



AF R 



IJ* 



(3.20) 



(3.21) 



where 

x d = (c a ,-s a ,isp,-icp), x s u = (s a ,c a ,-icp,-isp) 

for S° = if , h°, A and G°, respectively. In obtaining (3.20) we have set A d Y d to zero as 
done in refs. [27] and [28]. We will proceed in the same manner in deriving the formulae 
for the neutral Higgs couplings below. However, at the end of this subsection we will list 
for completeness the corrections to these formulae due to a non-vanishing A d Y d . 
It is now easy to check that the diagonal part (2.32) of the dS°d vertex is 

\x s ] JJ ee ( nf . mdj . + AF? S ) 33 = , mdj ^ (x s d + ~ejx s u ) . 



1 + e dj tan j3 



v d (l + ejtan/5) 



(3.22) 



For the flavour changing couplings 
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with J > I, neglecting S^ L and £y R and 



ignoring the contribution of A d Y d one gets from (2.34) 
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1 + e j tan (3~ u 
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^YVt^o 1 (x s u -x s d ttmp) 



, ! .7/ - -(I . •• (/ (3.23) 

v d (l + ejtanfj) 2 v ' 

in agreement with equation (9) of ref. [27]. Notice that effects of the genuine vertex cor- 
rection AFl s (proportional to ) are subleading compared to the effects of the rotation 
AD L oc Am d . In the case of split mass parameters TOq, the factor e Y in (3.23) depends 
on rrig of the J-th generation and mfj of the third generation. 
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Similarly, neglecting A^Y^ we find from (2.35) an explicit expression for 
J > I 
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eyjW(a;f -xf* tan/5) . 



(3.24) 



that has not been given in [27]. In particular this formula does not follow from the 



equation (9) of this paper, that for J > I is valid only for the 



X RL 



couplings. The 



expression (3.24) has been derived under the assumption that 6y(Qj,U^) ~ 6y(Qi,U[ 



Note, that in this case 



X LR 



ji 



can be obtained from 



x rl 



ji 



by simply replacing rridj 



by m dl and complex conjugation of x'^ and x a without any other changes in the indices. 
The couplings 



x rl 



ji 



and 



x lr 



ji 



for J < I can be found by applying the rule (2.36) 
to (3.24) and (3.23), respectively. 

In the same approximation it is easy to derive the relation between the CKM matrix 
Vji of the MSSM and the effective matrix Vjf. Using the hierarchy of quark masses and 
CKM matrix entries (and neglecting AU L ) one obtains from eqs. (2.21), (2.13), (3.12) 
the relations [24, 20, 27] 
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(3.25) 



We can now summarize the results for the neutral Higgs-quark vertices. In the case 
of .B-physics the pairs (J, I) = (3, 2), (3, 1), (2, 3) and (1,3) matter. Combining (3.23), 
(3.24), the corresponding formulae for J < J, and (3.25) and neglecting the dependence 
of e and ej on the generation indices we find the final formulae for the neutral Higgs 
couplings relevant for 5-physics that exhibit all large tan (3 enhanced corrections: 
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that are valid for both J > I and J < I. 

In the case of if-physics the pairs (J, I) = (2, 1) and (1, 2) matter and we find 
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It is evident from (3.21) and (3.26)-(3.29) that the flavour violating couplings of G° vanish 
in this limit. 

Finally, we would like to list the corrections to (3.20), (3.22), (3.23) and (3.24) that 
result from 

(A d Y d ) JI = -y dj (e;5 JI + eyy^ 1 ) (3.30) 
with eo, ey and ej replacing eo, ey and ij in A u Y d , respectively. We find: 

5 (AFf ) U = (ej6» + e Y yW) x s d (3.31) 

and 

(x s d + ejxl) -> (4 - ejx s d + ijx*) (3.32) 

in (3.22). In [X RL ] JI and [X'[ R ] JI the correction A d Y d produces the following subleading 
terms: 

6 1 X Rl} JI = n r~ d \ m2 ^o J t-ey + (e e Y - e Y e ) tan/5] x s d (3.33) 
v d [l + ej tan fjy 

S[X S lr} JI = - n^A^m^ o 1 hey + (e ey - e y e ) tan/?] *? (3.34) 

which are smaller than the x d tan/3 terms in eqs. (3.24), (3.23) but can be of the same 
order of magnitude as the x^ ones. The analysis of sec. 5 shows that the above terms are 
not numerically important. 



3.4 Charged Higgs-fermion vertices 



In the same approximation we can consider also the corrections to the charged Higgs 
boson and charged Goldstone boson couplings. In order to find explicit expressions for 
the charged scalar (H^) and charged Goldstone boson (G^) couplings defined in (2.39) and 
(2.40), we combine the effective Lagrangians (3.3) and (3.14). Performing the rotations 
(2.6) and (3.17) we find 
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Expressing H^* and H[ M> in terms of H + and G + as follows 
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we calculate, the couplings 
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in (2.39) and (2.40) and 



the parameters 6jj L , e^f 1 , e G f and e G f defined by eqs. (2.48) and (2.49). 

For G + couplings, using (3.36) and comparing terms in square brackets of eq. (3.35) 
with eqs. (3.9), (3.10), (3.16), (3.18) we find 
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Subsequently using the relations (2.7), (2.17), we easily find(y / 2/f = #2/ 'y/2Mw) 
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This compared with (2.49) implies 
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The vanishing of e G \ and e G f is easy to understand. Indeed, in the effective Lagrangian 
constructed before the electroweak symmetry breaking and defined by eqs. (3.3) and (3.14) 
the Higgs mechanism must operate in the standard way (e.g. in the gauge the tree 
level exchange of G + must cancel the gauge dependence of the amplitude generated by 
the W + exchange) which requires that the G + couplings to quarks be precisely as in the 
SM with the mass parameters corresponding to the tree-level quark masses in the effective 
theory (that is m d (u)) and that the CKM matrix in the G + couplings is the same as the 
one in the W + couplings (that is V ). The same argument shows that ejf and e G f must 
vanish in the SU(2) x U(l) symmetry limit also without the approximation of vanishing 
electroweak gauge couplings. 

For H + couplings, from (3.35) and (3.36) we find 
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(3.41) 



We observe that in the first line of eq. (3.41) A^Y M cannot be neglected as it is multiplied 
by tan/3. On the other hand A U Y U can be neglected and Ul,r can be approximated by 
unit matrices. Consequently we find for the u R H + dL vertex 
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Approximating next m Uj by m Uj and expressing V through V eG by means of (3.25), we 
find 

31 = W^w C ° S ^ I 1 " tan/9 ^ + e ' Yyl5 ^ " A ^ ' (3 ' 43) 
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where 
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2 2 eye'y tan 2 j3 
1 + e tan p 



+ 1 (J,/) = (1,3), (2, 3) 
-1 (J,/) = (3,1), (3, 2) 
otherwise 



(3.44) 



(3.45) 



This agrees with ref. [28], where the presence of Aj/ has been pointed out. 

In calculating the HJ H + d R vertex in (3.41) we observe that the effects of A d Y d are 
not tan f3 enhanced with respect to the tree level coupling and can therefore be neglected. 
The second term involving A u Y d is suppressed for large tan/5 with respect to the first 
term and can also be neglected. Consequently 
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Using the formulae (2.15), (2.16), (3.12), (3.25) and retaining only terms dominant for 
large tan/3 (and unsupressed by quark mass ratios or small CKM matrix elements), we 
find 
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= tan/3 V 

V2M W cff 1 + ejtan/5' " J1 l + ejtan/3 

for all (J, I) with ij defined in (3.6). 

Inspecting (3.47) we find that in the case of the vertices involving V^ ff and V t e d the 

corrections depend on the top Yukawa coupling yf while the corrections to vertices 

and do not. 

We would like to emphasize that whereas the rule (3.47) for (J ^ 3, 1) and accidentally 
for J = I = 3 is equivalent to expressing in the tree level formulae Vji and m dl through 
Vjf and fn dl respectively, for J = 3 and I — 1, 2 it is more involved. Expressing in these 
cases only Vji and m dl through Vjf and m dl) would give (1 + e 3 tan /3)/(l + e tan (3) 2 as 
opposed to the correct result 1/(1 + e 3 tan/3). 



3.5 Chargino-Quark-Squark Vertices 

As argued in subsection 2.6, to account for rotations (2.8) performed on the quark fields 
in the higher dimensional operators resulting from integrating out charginos, their Wilson 
coefficients should be calculated by using the modified chargino-quark-squark couplings. 
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Modifications of c l L given in (2.53) are trivial. In the SU(2) x U(l) symmetry limit 
the modifications (2.54) of c l R can be easily given in closed form if we note that they are 
exactly as in the expression (3.46) for the couplings [Plr] ji ■ Thus, in c R in (2.52) one 
should replace as in (3.47) 



m dl 



1 + e j tan p 

3.6 Introducing Flavour Dependence 



(3.48) 



Let us next generalize the main results of this section to the case of flavour dependent eo, 
ey, e' and e' Y - In this case (3.5) and (3.12) are replaced by 

(A u Y d ) JI = -y dj (e&V + eMytV'V* 1 ) = -y dj (lj5 JI + e^y^ 1 ) (3.49) 



(Am/ 1 = m dj tan/? (e,j5 JI + e? J) y t 2 A J/ 



where now 



e j — e + e Y y t o 



(3.50) 
(3.51) 



and it is ij given in (3.51) that enters the relation (3.13). 

In agreement with the typical pattern of the soft SUSY breaking mass parameters, in 
what follows we will assume that 
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With this assumption, instead of the relation (3.25) for (J I) = (13), (23), (31) and (32) 
we find: 
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with ej given in (3.51) and e defined through 
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Note that e differs from e J ^ that enters (3.51). 



Similarly we find for J > I 
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where 
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ej + {h - ej)(ey 7) / e r J) ) tan /^ 
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(3.57) 
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and 
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JI 



for J < I can be obtained by using the rules in (2.36). These 
improvements can be easily incorporated in the formulae (3.26)-(3.29) taking into account 
that in expressing A- 7 through Vj/, the formula (3.53) instead of (3.25) should be used. 
Next for the right-handed charged Higgs couplings we find 
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with e defined in (3.54). Finally 
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Flavour dependence in the formulae (3.43)-(3.45) for the couplings 
introduced in sec. 5.3. 
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(3.61) 
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3.7 Summary 



In this section we have derived explicit expressions for the flavour changing neutral 



scalar couplings 
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and 
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and 
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and the charged scalar couplings 

at large tan/3 in the SU(2) x (7(1) symmetry limit assuming the 
dominance of a s and the top and bottom Yukawa couplings. The main results are given 
in (3.23), (3.24), (3.26)-(3.29), (3.39), (3.40), (3.43) and (3.47). The important formula 
[20, 27] that relates V n to Vjf is given in (3.25). We have also indicated how these 
formulae should be changed to account for flavour violating effects in the case of split 
squark mass parameters. This set of formulae allow to calculate scalar contributions to 
any flavour violating process. In addition we have presented the rules for incorporating 
large tan (3 effects in the Wilson coefficients of the higher dimension operators generated 
in the course of integrating out sparticles. 
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4 Comparison with the Literature 



The formulae presented in sections 2 and 3 generalize and in certain cases correct those 
present in the literature. Let us discuss this in details. 

The results in section 2 are new as they go beyond the limit of zero electroweak gauge 
couplings and also beyond the SU(2) x U(l) symmetry limit. We devote the next section 
to the detailed comparison of the approximation of sec. 3 and the complete approach of 
sec. 2. Our results of section 3 agree with the recent analysis in [28] if one makes the 
following replacements in our formulae 



-e > e YVt -»■ e i, tyVb 



(4.1) 



with e unchanged. We would like also to compare our results with those of [25] in order 
to have a dictionary between the two notations. This will prove helpful for discussing 
B — > X s 7 decay. 



The explicit expressions for 



,H(G) 



.ii 



,ff(G) 



JI 



with 



RL with J = 3, I — 1, 2, 3 and for ± LR 
J = 1,2,3 and 1 = 3 given in eq. (17) of ref. [25] do not take into account the modifications 
of the CKM factors summarized in (3.25). However, correcting the formula (17) of [25] 
has to be done with care. In particular, identifying V used in that paper with V c q would 
generally give wrong results for 



D H(G) 
LR 



JI 



while it is correct for 



D H(G) 
RL 



JI 



if the term 



Ajj given in eq. (3.45) is neglected. 

In order to see this explicitly let us observe that the couplings 



D H{G) 
RL 



JI 



involve left- 



handed down quarks that enter also the W -quark couplings. As the corrections resulting 

after 



D H(G) 
RL 



JI 



from rotating the ul and ur fields can be neglected, the CKM matrix in 
the rotations (2.8) becomes automatically V eS except for V in the last term in eq. (3.42). 
Approximating V in the latter term with V we find P RL in (3.43) with Ajj set to 

zero. This result, obtained in the SU (2) x U(l) symmetry limit, agrees then with P RL 
in equation (17) of [25] provided the CKM matrix there is identified with V cS . We find 
also the relation between ej/ introduced by us and the parameters e'j(I) in [25]: 



,HL 
Z JI 



e'j(I) tan/3 



J = t, I = d,s,b, (Aji^O) 



(4.2) 



However, as emphasized in ref. [28], Ajj cannot be generally neglected for |eytan/3| and 
\e' Y tan(3\ larger than 0.5 and it could be important for e' ks —e' Y when the C(tan/3) term 



in ej/ is small. 



However, we will see in section 5 that even with the Aji term included, €jj L obtained 
in the approach of sec. 3 deviate rather significantly from ej/ computed as in sec. 2. The 
replacement of e' + e' Y yl5 n in eq. (3.44) obtained in the strict 577(2) x f/(l) symmetry 
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limit and neglecting couplings other that a s and Yukawas by e'j(I) given in eq. (16) of 
[25] does not help much. To get a reliable approximation one has to include additional 
terms that depend on the electroweak gauge couplings which are numerically important 

r „ -I J I 

but have been neglected in all previous analyses. In contrast to P^ L , the coupling 

Pr L of the charged Goldstone boson as used in ref. [25] agrees with our if V of [25] 
is interpreted as V e Q. As we will see, e'jj computed in the approach of sec. 2 are indeed 
negligible. 



In the case of 
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JI 



the right-handed down quarks are rotated differently from the 



left-handed down quarks present in the W^-quark couplings and this mismatch results 

r h(g)~\ 

in additional tan/3 enhanced corrections. Interpreting V in \Pj^ji in equation (17) 
in [25] as V eS , would generally miss them. However, it turns out that for J = 1, 2, 3 
and 1 = 3 considered in [25], the results correct in the SU(2) x U(l) symmetry limit 
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of this paper not with V e s as in the case 



but with the original CKM matrix V. Indeed, expressing the V elements 



ji 



of [25] in terms of V by using the rules in (3.25), the formulae for these 
couplings given in [25] can be transformed into ours provided the following replacements 
are simultaneously made: 
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in the case of 
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and 



V JI l + e' 7 (J)tan/3 
Vjf l + e/tan/3 
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in the case of 
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(note that ej of [25] corresponds to ej in our paper). As in the 



SU(2) x U(l) symmetry limit 

e' b (t) 



£b — £3, 



e 0, 



(4.5) 



with €3 and eo defined in (3.6), the formulae given above reduce indeed to the correspond- 
ing formulae for the parameters €jj G ^ R given in (3.40) and (3.47). It should be emphasized 
that the factor V JI /Vjf in (4.4) that we introduced here is essential for vanishing of e^p 
in the SU(2) x U(l) symmetry limit as required by general arguments. We will see, how- 
ever, in the next section that for split squark masses the formula (3.47) for e^/ 2 has to 
be modified along the lines of section 3.6 in order to approximate well the results of the 
complete calculation. We shall also explain in more detail the status of the recipe (4.4). 
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5 The Size of the Large tan (3 Corrections: Full Cal- 
culation and the SU(2) X 17(1) Limit 

This section is devoted to the numerical evaluation and the assessment of the importance 
of the corrections discussed in sections 2 and 3. We will also compare results of the 
complete diagrammatic calculation of these corrections with the approximation based on 
the 577(2) x U(l) symmetry limit discussed in section 3. This will allow to test the quality 
of the latter approach and to find its limitations. We will also present analytic expressions 
for various couplings that approximate the full approach much better than the formulae 
of section 3. 

5.1 Full Calculation 

The one-loop contributions of sparticles to quark self-energies S^ L( -^(0), £y L( ^(0) and 
various vertex formfactors AFl(r) at vanishing external momenta, that are the main 
ingredients of the full method of section 2, can be easily calculated by using the Feynman 
rules of the MSSM collected in [30] and the standard one-loop functions. We recall the 
relevant expressions in the Appendix A.l. 

In order to find the original mass parameters rridj and the original CKM matrix V JI 
that enter these rules we proceed as follows. We begin by computing defined in (2.30). 
Since to a good approximation (S^ L ) JJ oc rridj the quark mass cancels out in the ratio in 
(2.30). Therefore, in the actual calculation of the MSSM mass parameters m^ 7 can be 
replaced by rridj that are known experimentally. Note also, that to compute diagonal self 
energies one can use experimentally known V oS instead of the yet unknown V. Having 
obtained e dj in this manner we can calculate rridj ~ rridj/ {1 + tan/3). 

In order to find V we proceed iteratively: starting with rridj and V = V we compute 
A ; m d by means of (2.5) and (2.11) and subsequently using (2.13) and (2.14) we find 
the rotations T>l,r = 1 + AD^. This allows to find V in the first approximation by 
inverting the relation (2.21) in which AUl can be set to zero. This is next used to 
compute A'm d anew and so on. The procedure converges quickly yielding V and the 
rotations D^. Having these, we can compute all the necessary formfactors AFl,r and 
the couplings X RL , X LR , P RL and P LR . 

The simple and transparent expressions derived in sees. 3.2-3.5 are very convenient, 
as they make the effects of large tan/3 explicit and allow for easy qualitative discussion 
of the dominant effects. However, as the comparison with the complete calculation will 
reveal, they are not very accurate numerically and it is desirable to improve the accuracy 
of the approach of section 3 retaining, however, its simplicity. It turns out that for the 
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CKM matrix and the couplings ^Crl], [-Xl.r], and [Plr] this can be achieved by simply 
using the expressions (3.53)-(3.61) with tj and ey 1 ^ calculated directly from the full one 
loop correction Am^ to the down-quark mass matrix given in (2.5): 



^) JI for J I, ~ej - {Amd)JJ . (5.r 
rridj Xq 1 tan (3 rridj tan (3 



The complete formula for (Am c i)' JI is given in (2.5) with the relevant expressions for 
quark self-energies given in the Appendix A. 2. Note that with Sy ^ and ij all parameters 
entering (3.53)-(3.61), in particular e in (3.54), can be calculated. In this manner also 
the effects of the electroweak couplings g\ and g 2 present in (Am d ) JI are automatically 
taken into account. Simple but accurate expressions for the couplings [Pr L ] can also be 
given. They will be discussed separately below. 

5.2 The neutral Higgs boson sector 

In order to implement the formulae of section 3 we need the quantities eo, ey defined by 
eq. (3.5) and e' , e' Y of eq. (3.15). Evaluating diagrams shown in fig. 5 we find 

where x Q/g = mq/m-, x D/9 = m 2 D /mj, x®^ = m 2 Q /\^,\ 2 etc., and rag, m 2 D , mjj, A t , and 
A b are the parameters of the soft supersymmetry breaking in the MSSM Lagrangian. 3 
The function H 2 (x,y) is defined as 

u( \ xlnx i ylny (KA\ 

H * X < v)= il- X )i X - V ) + (l-y)<y- X )- <5 ' 4) 

As discussed in sec. 3, for squark mass parameters itlq, m 2 D , mfj having different values 
for different generations (e.g. if the the third generation squarks are split from the others 
as happens in the minimal SUGRA scenario) eo, ey, e and e' Y depend on the generation 
indices. We have given the relevant expressions in subsection 3.6. In our comparison 
of the approximate and full methods, we parametrize the squark spectrum with four 
independent mass parameters: m 2 = 2 = m 2 Dl 2 = m 2 j i 2 for the squarks of the first 
two generations and 771q 3 , m 2 D3 , m^ 3 for the third one. Consequently, in eqs. (3.26)-(3.29) 

3 Our convention for A u and Ad parameters is fixed by the form of the left-right squark mixing terms 
in the squark mass squared matrices which read — m u (A u + fi cot 0) and — md{Ad + ^tan/3) for the up 
and down squarks, respectively. 
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€3 always depends on the third generation squark masses but as in sec. 3.6 we distinguish 
6q arising from expressing Vjj in terms of Vjf , which also depends on the third generation 
squark masses, and eo that result from expressing m dl2 in terms of m dl 2 which depends 
on m 2 . 

The formulae (5.2) for eo and ey apply in the strict limit of SU(2) x U(l) symme- 
try supplemented by the approximation g 1 = g 2 = 0. Under these conditions only the 
higgsino of mass couples to quarks through Yukawa couplings. In application to the 
relation between m dj and m dj and to the neutral Higgs boson couplings it is easy to take 
into account the contribution of the true charginos. Simplifying appropriately the full 
expression for {T 1 f nL ) JJ (or the full expression for the one loop vertex diagram in which 
H° or h° couples to the up-type squarks and the charginos connect the external d-quarks) 
and comparing with the ey contribution to Am d (to AFl s ) one gets 

where xf^ C = mg/m^ etc., and the matrices Z + and Z_ are defined in ref. [30]. Similarly 
one could take into account also the left-right mixing of squarks. The formula (5.5) goes 
beyond the 577(2) x U(l) symmetry limit in that it takes into account the mixing of 
higgsinos and gauginos which depends on vacuum expectation values and on the gauge 
couplings g\ and g 2 and can be suitably extended to include also the mixing of left and 
right-handed squarks proportional to VEVs. However, it still neglects g\ and g 2 in the 
couplings of charginos to quarks and the contribution of the neutralinos. 

The magnitude of eo and ey can be easily estimated from eq. (5.2) remembering that 
H 2 (l, 1) = —1/2 and that H 2 (x,x) varies between —0.18 for x = 4 and —1.7 for x = 0.1. 
Figures 6a,b show the comparison of e d2 ~ e dl and of e d , A computed from the full one-loop 
self energies (eq. (2.30)) with eo and e^, respectively, for a sample of points in the MSSM 
parameter space corresponding to sparticles heavier than 500 GeV. As seen in figs. 6a 
and 6b, typically |e<2 3 | and \ed 2 \ (|eo| and | ^3 1 ) are of order ~ 5x 10~ 3 and reach ~ lx 10~ 2 
only for very special combinations of squark and gluino masses, and/or \A t \. 

The two distinct bands seen in panels a) and b) of the figure 6 correspond to // < 
(e < e d2 < and e 3 < e^ 3 ) and fi > (e > e d2 > and e 3 > e d3 ). The discrepancy 
between eo (is) and e d2 (e d3 ) (quantified by the deviation of the bands from the diagonal) is 
partly due to the absence in eo given by (5.2) of the neutralino contribution to the d-quark 
self energies but mostly to the fact that in the full approach the chargino couplings depend 
also on the gauge coupling constants so that charginos do contribute also to e and hence 
also to e dl and e d , 2 . Neglecting the contribution of A d Y d in the formula (3.10) corresponds 
in the full approach to neglecting the contribution of A b (of ji) to the mixing of the left 
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Figure 6: Comparison of the results obtained in the complete calculation with the ones 
obtained in the approximation based on the SU(2) x U(l) symmetry limit for a sample 
of the sparticle mass spectra for tan f3 — 50 and M A = 200 GeV. See the text for further 
explanations. 
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and right sbottoms (stops) and is justified for \Ab\ <C tan/5 and \fi\ <C A t tanf3. More 
detailed comparison shows that the differences between eo (£3) and e^ 2 (e^) are of order 
15%-20% when the epsilons assume phenomenologically relevant values. 

The gap in values of €d 2 (eo) seen in panel a) of fig. 6 reflects the existence of a lower 
bound on [mig in the scan which follows from the imposed condition mc, , m g > 500 GeV. 
Note also, that while e and e Y as given by (5.2) or (5.5) are UV finite, have to be 
renormalized in the MS scheme. We have checked however, that the variation of e dl with 
the renormalization scale Q is very small for 500 GeV< Q < m g . 

Figure 6c shows the comparison of 4 V t °f /V td obtained from the complete expression 
(2.21) and in the approximation (3.25) for a sample of sparticle parameters and tan (3 = 50. 
It is to be noted that the approximation (3.25) with e and ty computed as in eqs. (5.2) 
and (5.5), works only qualitatively. 

Figures 6a - 6c show that in general the effects of the electroweak couplings are non- 
negligible. They could be included in the framework of the SU(2) x U(l) symmetry 
limit by computing diagrams with gaugino exchange but since the full expression for the 
supersymmetric contribution to the down quark self energy are simple enough, it is more 
practical to use the formulation of section 3.6 with ey ^ and e j computed as in (5.1). The 
formula (3.53) improved in this way gives perfect approximation to the ratios V t c f /V t d etc. 
obtained in the approach of sec. 2. Obviously, ij given by (5.1) coincide with the factors 
Ed j of sec. 2. 

The panel d) of figure 6 shows for Ma = 200 GeV and tan/5 = 50 the correlation 5 of 
BR(B — > Xs'-f) with the factor (1 + eotan/3)(l + 63 tan (3) that enters the denominator 
of the flavour violating neutral Higgs boson couplings (3.26) and (3.27) improved as 
in sec. 3.6. We have used here the eq. (5.1) to determine eo and e 3 . We observe 
that the experimentally acceptable value of BR(B — > Xs'-f) can be realized both for 
(1+eo tan (3) (1+63 tan f3) greater than one and smaller than one corresponding respectively 
to the suppression and the enhancement of the actual scalar FCNC compared to the naive 
one loop estimates of ref. [6, 22, 23]. 

In contrast, it is interesting to observe in the panel e) of figure 6 that the correlation 
between V t e d s jVtd and B — > X s "y has quite different implications. Here the requirement of 
acceptable rate of the B — > X s ^ decay eliminates (for Ma = 200 GeV and tan (3 = 50) all 
points with V t f /V td < 1. 

Both facts can be qualitatively understood by investigating signs of e and e^ 1 -* and 
remembering that, as follows from eqs. (5.2), e ~ e oc \i and e Y ~ ey oc — fiA t 
4 We have checked that the full approach based on eq. (2.21) also gives Vff /V t d ~ V t c s s /V ts ~ 

vtf/v* « v c f/v cb . 

5 For details of the BR(B — ► X s -f) calculation see sec. 6. 
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(the function H2 is negative). Let us consider the case eo < 0, £3 > which according 
to eq. (3.53) leads to Vf d < jVu < 1- This requires fi < and A t > 0. However, the 
relative sign of the chargino-stop and H + contributions to the amplitude of the B — > X s 7 
decay is opposite to the sign of jiA t . Therefore, jiA t > is a necessary condition for the 
cancellation to occur [26]. Hence, for e < 0, e 3 > no cancellation between the chargino- 
stop and H + contributions is possible which explains why V t e f /V t d < 1 is incompatible 
with the BR(B — > X s ^) constraint. 

On the other hand, for (1 + e tan/3)(l + e 3 tan/3) < 1 (enhancement of the flavour 
changing scalar couplings) that holds always for fi < such a cancellation is possible 
for A t < 0. This is further illustrated in figure 6f which shows the correlation of e 
and e 3 for those sparticle parameters which for tan (3 = 50 and M A = 200 GeV give 
2.81 x 10~ 4 < BR(B -> X s 7) < 3.65 x 10~ 4 . It should be stressed that the frequently 
used argument [27] that obtaining acceptable B — > X s ^ rate requires /i > 0, which implies 
the suppression of the scalar FCNC, applies only in particular SUSY breaking scenarios. 
For example, in the minimal SUGRA, due to the fixed point structure of the relevant 
RGEs, not too large values of A (common value of A t , A b etc. at the high scale) always 
lead to a fixed sign (positive in our convention) of A t . However, fixed sign of A t needs 
not be a common feature of all SUSY breaking scenarios and can be reversed even in 
the minimal SUGRA if the initial value \A \ is sufficiently large (see e.g. ref. [32]). It 
is therefore justified to study the phenomenological consequences of both possibilities: 
(1 + eotan/3)(l + e 3 tan/3) greater than one and smaller than one. 

Finally, let us stress that in the case of the enhancement of the scalar flavour changing 
couplings the value of (1 + e tan/3)(l + e 3 tan/3) with e and e 3 calculated from eqs. 
(5.2), (5.5) is smaller than (1 + eo tan/3)(l + e 3 tan/5) with e and e 3 calculated using (5.1) 
by 10%— 35% which means that the naive approximation gives too big an enhancement. 
Similarly, in the case of the suppression of the scalar flavour changing couplings the value 
of (1 + eotan/3)(l + e 3 tan,3) obtained by means of (5.2) and (5.5) is larger by 10%— 25% 
than the corresponding factor calculated using (5.1). Thus the approximation of sec 3. 
overestimates the effects of resummed large tan f3 corrections in the neutral Higgs boson 
couplings. 

This is reflected in figure 7a where we show the comparison of the couplings 
H° to the down quarks obtained from the complete calculation and from the approximate 
formulae (3.26) with e and ty computed by means of eqs. (5.2), (5.5). In this figure 
positive (negative) values of X^ L correspond to negative (positive) value of the product 
fiA t . Clearly visible shorter band corresponds to positive /i (suppression of the flavour 
changing scalar couplings by the factor (1 + e tan/5)(l + e 3 tan/3)) and the longer one to 

r 1 32 

ji < (enhancement). For points for which the coupling X^ L has a non-negligible value 
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Figure 7: As in figure 6 but for the couplings [X RL ] JI of H° to the down quarks. 



the difference between the full and approximate calculation varies between 10%— 25% for 
pi > (suppression) and 10%— 40% for /i < (enhancement). Other flavour changing 



couplings: 



Y H 



n31 



X RL 



1 23 



etc. exhibit similar behaviour. For completness, in figure 
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X RL 



of H° obtained from the 



7b we show the comparison of the diagonal couplings 
complete calculation and from the approximate formulae (3.22). 

Similarly as in the case of the ratio V{f /V td) the quality of the approximation of section 
3 seen in figure 7 can be greatly improved by using the formulae (3.22) and (3.26), (3.27) 
improved according to the rules of sec. 3.6 with e and ij and Sy ^ computed as in eqs. 
(5.1). Since the tan/3 enhanced parts of the flavour changing couplings are proportional 
to A t , such an approximation is not good only if the parameter A t is small but then also 



the couplings 



Y s 



JI 



Y s 



JI 



for J 7^ I become phenomenologically uninteresting. 



For (almost) degenerate squarks of the first two generations the approximation of 



subsection 3.6 for couplings 



A RL 



and 



Y s 
^LR 



LI 



where 1 = 3 and/or J = 3 with e and 



I j and €y ^ given by eqs. (5.1) could break down also for unrealistically large differences 
m 2 — itiq 3 . This is because, as we have explained in section 3.1, we are working in the 
scenario in which the left-left block of the up-type squark mass squared matrix is VrnqV^ 
and is therefore not diagonal in the generation space as we allow for mg 3 ^ m 2 = = 
mQ 2 . However, due to the special structure of the CKM matrix this flavour violation 
remains negligible as long as \(nT,Q 3 — m 2 )V t bV* b \ remains small compared to mq 3 . If this 
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condition is not satisfied then, the couplings of stops to charginos and down quarks depart 
from the simple structure (2.52) and in consequence (A'm ( i) JI is no longer proportional 
to Aq 7 . Since in our scans \(mQ 3 — m 2 )V t bV* b \ <C mL the effects of this additional flavour 
violation never manifest themselves. The main effects of flavour violation arise then from 
flavour dependence of e and ey and when the approximate formulae (3.26)-(3.29) are 
modified according to the prescriptions of section 3.6 with epsilon parameters calculated 
by means of (5.1) they always agree very well with the results of the full calculation. 

5.3 Charged Higgs Couplings 

In this subsection we will asses the magnitude of the factors e^/, Cjf" and ej/\ e G jf intro- 
duced in subsection 2.5 to parametrize the corrections to the charged Higgs and Goldstone 
boson couplings and test the quality of the approximations developed in subsection 3.3. 
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Figure 8: Dominant corrections to the charged Higgs boson couplings to quarks. 

The factors e' Q , e' Y , given in eq. (5.3) entering the expression (3.44) for ejj L obtained in 
the SU(2) x £7(1) symmetry limit and for vanishing electroweak gauge couplings can be 
easily extended beyond this approximation. Improvement analogous to the replacement 
of ey of eq. (5.2) by ey of eq. (5.5) in the neutral Higgs boson sector has been given in 
ref. [25]. Simplifying appropriately the full MSSM expressions for the diagrams shown in 
fig. 8 allows to replace [i by the contribution of the four neutralino mass eigenstates and 
to take into account the mixing of the left and right squarks of the third generation. As 
a result e' + e' Y yl in the formulae (3.15) and (3.44) for (A d Y u ) JI and ejf, respectively 
are replaced by the flavour dependent quantity e'j(J) of [25]: 
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where x^ 9 = M| /m| (xt /9 = Ml Am§), k = 1,2 with w fc (c4) being the up (down) 
squark of the J-th (7-th) generation and and cj etc. denoting cosines and sines of the 
appropriate mixing angles. (For the first two generations c-j ~ 1, s- j ~ 0.) Similarly, 



„<V<? 



a; 



u/AT 



= M| Jm 2 N , = 1, 2, Z = 1, 4 etc. In the formula (5.6) ?/& is given by 



y/2 



(5.7) 



Using the identity Z^m N[ Z^ 
e'j(I) ^ e' + e' Y yl 



v d (1 + e 3 tan/7) 
-/i* it is easy to see that in the SU{2) x {7(1) limit 




Figure 9: Additional important corrections to the charged Higgs boson couplings to 
quarks. 

Similarly as in the case of the replacement of ey of eq. (5.2) by that of eq. (5.5), the 
use of e'j(I) does not improve the approximation to €jj L obtained in the full approach 
of sec. 2.5. This is because while going beyond the ST/ (2) x {7(1) symmetry limit in 
taking into account vacuum expectation value dependent sparticle mixing, it still neglects 
contributions which are present in the 577(2) x £7(1) symmetry limit for nonzero gauge 
couplings g\ and g 2 . The latter effects can be incorporated taking into account also the 
SU(2) x i7(l) symmetry breaking effects by computing diagrams shown in fig. 9 (in the 
SU(2) x {7(1) symmetry limit they correspond to diagrams shown in 10 contributing to 
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ArfY„ for nonzero g± and g-z). Taking into account their dominant parts (which are finite) 
amounts to adding to e'j(I) the following terms: 



x J m - 1 & yyH^iJ^ff, ( r N / c r Q/°) 
6atj[I) ~ 16** ^2c w \ Z + ZNa m Cl H2 <^ ' Xjl ) 

where = m 2 N ./m 2 Ci) xj( C = rnQj/rriQ, etc. and a'- 7 = Z^(swZ^ + cwZ 2 ^) 



I I 




a) b) 

Figure 10: Additional contribution to A^Y„ for nonzero electroweak gauge couplings. \ 
and if) denote U(l) and SU{2) gauginos, respectively. 

Only the sum of all these corrections including the delta term (3.45) reproduces the 
correction €jf computed as in sec 2.5: 

e H jt « tan/3 [e'j(I) + 5 a e'j(I) + 5 b e'j(I)} + A n (5.8) 

While in most cases the term proportional to a s in e'j(I) is dominant, the term 5 a e'j(I) + 
5b€j(I) is usually more important than the term Ajj and even than y\ dependent part 
of the e'j(I) term. Let us also stress that the quality of the approximation of ejf by the 
rhs depends crucially on how well €3 approximates the factor e<f 3 (2.30) needed to obtain 
yb and how well epsilons reproduce the ratio V JI /V^ . Since the approximations of sees. 
3.2-3.3 are usually not better than 15%-20%, the contributions to the lhs and rhs of eq. 
(5.8) not proportional to a s usually differ by about 20%-50%. Reliable approximation 
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can be obtained only by using in (3.45) eo and e Y computed as in eq. (5.1) instead 
of eo and ey. The approximation (5.8) works then to an accuracy 5%-10%. This is 
demonstrated in figure 11 where in panels a and b (panels c and d) we compare e 33 L 
( e 32 L ) calculated as in section 2.5 with the approximation based on the formula (3.44) 
and with the approximation based on eq. (5.8), respectively. In panel a (and c) we have 
used in the Ajj term (3.45) eo and ey given by eqs. (5.2), (5.5) whereas in panel b (and 
d) e and ey 2 - 1 from the formula (5.1). The difference between the results of the full and 
approximate calculations is quantified by the deviation of the points from the diagonal 
line. The dramatic improvement of the approximation is clearly seen in panels b and d 
of figure 11. 




Figure 11: As in figure 6 but for the corrections e 33 L and e 32 L to the couplings of H + to 
the left-handed down quarks. Only points yielding acceptable BR(B — > X s 7) are shown. 

Next we consider the correction e^j R . As in the approximation based on the SU(2) x 
£7(1) symmetry limit, also in the complete calculation of sec. 2.5 it is entirely determined 
by the rotations in the down quark sector (the genuine vertex corrections and the con- 
tribution of the vector self energies are negligible) and can be obtained from the formula 
(3.46) with exact matrices and D^. The quality of the approximation (3.47) depends 
therefore crucially on how well the matrices D L and D R are approximated. For split 
squarks the full flavour dependence of the correction (Am d ) JI is due also to the flavour 
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dependence of eo and ey and, consequently, the matrices and D# obtained in the strict 
577(2) x U(l) symmetry limit with universal eo and ey deviate from the true ones. To 
obtain a reliable approximation for €jP in the more realistic situation, in which squarks 
from the third generation are split from the others, one has to use flavour dependent 
expressions (3.58)-(3.61) and the improved formulae in (5.1). The results of this improve- 
ment are shown in figure 12 where we compare e^ 2 R obtained in the full approach based 
on eqs. (2.42), (2.44) with the approximation (3.47) using universal e and ey obtained 
from eqs. (5.2) and (5.5) (panel a) and with the ef 2 R obtained from eqs. (3.61) and (5.1) 
(panel b). 
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Figure 12: As in figure 6 but for the corrections e^ to the couplings of H + to the 
right-handed down quarks. Only points yielding acceptable BR(B — > X s 7) are shown. 

Finally we consider the corrections of the charged Goldstone bosons to quarks defined 
in eq. (2.49). As we argued, in the SU{2) x U(l) symmetry limit e^f = e G 3 f = 
also for nonzero electroweak gauge couplings gi and g 2 . To understand our results, it 
is important to stress at this point that integrating out sparticles and constructing the 
effective Lagrangian in the 577(2) x U{1) symmetry limit without making any assumption 
about dominance of some particular couplings (like a s , or y t ) requires also the inclusion 
of the vector parts of the quark self energies and of all vertex diagrams, in particular the 
ones shown in fig. 10. In other words, the whole procedure of sec. 2 has to be used 
with the exception that the Higgs boson vacuum expectation values are equal zero. The 



41 



argument for vanishing ejf and e^P applies then up to 1/M| USY corrections. 

To go beyond the SU{2) x U(l) symmetry limit we write the factors ejf and €jP as 

JK 



V2M W [U l R ■ [AFt ■ VI - ^Am.) • U L j V« 



92 rn Uj V^ 



V2M W [Dt ■ (yt . A po + ^-(Am d )t) • V R 



e Ji - 

GR - n ifJI^r 



KI 



92 V&m dl 



(5.9) 



(5.10) 



where AF^ R are defined in eqs. (2.43), (2.44). (Similar expressions can be written 
down for corrections ejj L and ejj R .) As can be seen by using the formula (3.35), in 
the SU(2) x U(l) symmetry limit AFf = ^-Am u V, AF R = -^^(Am^ and 
the corrections indeed vanish. (It is also easy to see that Sy L and ^y R cancel out 
in (5.9) and in (5.10), respectively.) Because both, V^AF R and [g2/ \f 7 2My/)( y Am d )^ 
contain contributions oc tan/3, any mismatch between them could in principle give rise 
to tan/5 enhanced terms in ejP and this motivated the authors of ref. [25] to keep the 
correction e'jf in their BR(B — > X s ^i) calculation (no such enhancement was expected 
in €jj). Exact numerical evaluation of the formulae (5.9), (5.10) reveals, however, that 
the cancellation between AF^ ■ and —(g 2 /V2M w )Am d in (5.9) and between V* • 
AF R and (g 2 / V2M w )(Am d y in (5.10) is surprisingly accurate when these quantities are 
calculated as in sec. 2 which means that the £77(2) x U(l) symmetry breaking corrections 
to AF^ R and Am Utd are tiny compared to the corrections introduced by switching on the 
electroweak gauge couplings. We find that |ejf| and \ej^\ are of the same order of 
magnitude and are ~ O(10~ 4 ) in most of the parameter space (for sparticles heavier than 
500 GeV and tan/? = 50). These are small corrections compared to 1 and do not influence 
appreciably the phenomenology we are going to discuss. Note also, that approximating 
Am d or AF^ R by dropping some "small" terms or using simplified couplings in the vertices 
as compared to the ones given in ref. [30] spoils this delicate cancellation and generates 
\e^' R \ ~ lO-M ! 

Finally let us explain the status of the rules (4.3) and (4.4). As long as one considers 
only 0(a s ) corrections to the effective Lagrangian couplings the rules (4.3) and (4.4) 
are the true equalities. This has to be so, because when the charged Goldstone boson 
vertices corrected as in [25] are used in the one loop expressions for the Wilson coefficients 
CV 5 g of the (chromo) magnetic operators they reproduce correctly [26] appropriate tan/3 
enhanced terms in the two loop contributions to calculated in [33]. The ideology of 
the effective Lagrangian serves then only to correctly resumm these terms. Indeed, in this 
case V eS = V, Dl,r — I and inserting AF R in the form 

(AF£) JI = --j9_V JI m il 4(J)t f m0 (5.11) 



42 



as in ref. [25] in (5.10) one recovers the rule (4. 4). 6 However, when the corrections 
depending on the top Yukawa coupling to the vertices of the effective Lagrangian are taken 
into account, the true formula (5.10) for e^P is more involved and cannot be obtained 
from the rule (4.4) just by including the terms oc y\ to e' 7 (J) as in [25]. Formulae (4.3) 
and (4.4) are then only the substitution rules allowing to translate the formulae of [25] 
into ours. 



6 AM d , s , 



B d,s 



and B — > X s ~f 



The B% 



fi + fi~ decays and the B° s -B° s mixing attracted recently a renewed attention 
due to the observation that for large values of tan/3 their amplitudes can receive very 
large contributions from diagrams depicted in figures 13 and 14 in which the black blobs 



represent the flavour changing couplings 



,11 



and 



^LR 



.11 



discussed in the preced- 



ing sections. These contributions have been found to increase by orders of magnitude 
the branching ratios for the rare decays B® d — > [20, 6, 21, 22] and to decrease 

substantially the B° s -B° s mass difference AM S [23]. 




tan 2 8 # «( tan 8 

h°,H°,A° ^ 



Figure 13: Diagrams giving dominant contribution to B® d — > l + l amplitudes at large 
tan/3. 

As demonstrated in [20, 6, 21, 22], for tan/3 ~ 50 and non-negligible values of the 
parameter A t the B° sd — > fi + n~ amplitudes are totally dominated by the diagram of 
fig. 13. In the absence of any other constraints on the MSSM parameter space, the 
corresponding branching ratios, which behave as | A t tan 3 /3/M\\ 2 , can be enhanced by 
up to three orders of magnitude relative to the SM predictions and can even exceed the 

6 However, calculating e'j(J) from the triangle diagram of fig. la of [25] with simplified G + couplings 
to up and down squarks would result in e'jf ~ 10~ 2 , much bigger than O(10 -4 ) obtained with the full 
G + couplings given in [30]. In particular, we have found that for non-negligible sbottom mixing (if, say, 
\Ab\ ~ \At\) neglecting in this vertex terms oc y\ and cx g\ has dramatic effect on the cancellation in 
(5.10). On the other hand, neglecting simultaneously the contributions of the vector self energies to AF^ 
and to Arrid does not affect it. 
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present experimental bounds: 



BR(B° d -> /i+fj,-) < 2.1 x 1(T 7 BaBar [34], (6.1) 
BR(B° S -> < 2.0 x 1(T 6 CDF [35] . (6.2) 

The B a sd -B a sd mass difference in turn receives several contributions: 

AM S = |(AM S ) SM + (AM/" + (AM S ) X± + (AM S ) DP | = (AM S ) SM |1 + f s \ (6.3) 

(by definition AM S is a positive definite quantity). For large tan (3 and non-negligible 
A t the contribution (AM S ) DP of the double scalar penguin (DP) shown in fig. 14 is 
the dominant correction to the SM contribution (AM S ) SM . Both, (AM S ) DP and the 
contribution (AM s ) ff± of the box-diagrams with top and charged Higgs bosons H ± have 
the signs opposite to (AM S ) SM . While generally smaller than (AM S ) SM , their sum leads 
for large tan (3 to a significant decrease of the predicted AM, (i.e. to f s < 0) independently 
of the choice of supersymmetric parameters. We will asses the relative magnitudes of the 
double penguin, charged Higgs and chargino box diagrams in sees. 6.1.3, 6.1.4. 




b R s L b L s R b R s L 



Figure 14: Double penguin diagrams contributing to AM, 



As is evident from the comparison of figures 13 and 14 there must exist a strong 
correlation between the enhancement of BR(B®^ —> fi + n~) and the suppression of AM S . 
In particular for < (1 + f s ) < 1 the experimental lower bound AM S > 15/ps puts 
an upper bound on the possible enhancement of BR(B®^ — > /i + /i _ ). Of interest is also 
the case (1 + f s ) < corresponding to a very large negative (AM S ) DP (that can be 
realized for very special values of supersymmetric parameters) which has quite different 
implications than the case < (1 + f s ) < 1. The main result concerning this correlation 
has been presented in [29]. Our purpose now is to explore it in more detail investigating 
in particular its dependence on the MSSM parameters and elucidating the impact of 
the tan j3 enhanced corrections to AM S on the standard Unitarity Triangle (UT) analysis 
which is necessary to determine the CKM matrix element Vu needed to predict accurately 
BR{B»^^^). 
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Large tan/5 effects manifest themselves also in the transition B — > X s 7, which being 
already relatively well measured, plays an important role in constraining the allowed region 
in the space of supersymmetric parameters [25, 33, 26]. Exploring the correlation between 
the increase of BR(B° S ^ — > fi + fi~) and suppression of AM S it is therefore important to 
incorporate this constraint in the analysis. 7 

In this section we collect the formulae for BR(B° S ^ — > /i + /x~) and (AM S ) DP in the 
large tan (3 limit and discuss the importance of various contributions. We recall also the 
UT analysis and include a subsection devoted to the calculation of BR(B — > X s ^). As in 
[23, 29] we will concentrate on the scenario with heavy sparticles and the mass scale of 
the Higgs sector close to the electroweak scale but will comment also on possible effects 
of lighter sparticles. The numerical analysis of the correlation of the enhancement of 
BR(B®, d j — > /x + /x~) and the suppression of AM S will be the subject of the next section. 

6.1 Formulae for calculating AM S 
6.1.1 The Effective Hamiltonian 

The effective Hamiltonian for B Q S -B Q S mixing has the form as follows 

r 2 M 2 

n g =2 = b^tygyuy E . (6 . 4) 

" i 

Here Qi are the AB = 2 operators and Gp is the Fermi constant. It should be stressed 
that with the improvements in the calculation of Higgs boson vertices the factorization of 
the CKM elements in (6.4) is only approximate. Still, it is an excellent approximation to 
the full calculation and allows for a transparent discussion of the unitarity triangle (UT). 
The set of dimension six AB = 2 operators consists of 

Q VLL = (hl,s L ) (Tl^sl) 

Q\ K = (b^s L ) (hs R ) (6.5) 

Qi L L = (fc) (fe) 
Ql LL = (5*vs L ) (^ u s L ) 

where = \ [7^,71/] and the colour indices are contracted within the brackets. The 
additional operators Q VRR , Qf RR and Q| RR are obtained from Q YLL , Qf LL and Ql LL by 

7 Large tan/3 effects have been afso investigated in (g — 2)^. However, they depend on the slepton 
sector parameters which can be correlated with the squark sector ones only within a particular scenario 
like e.g. minimal supergravity [36] or gauge mediation. 
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replacing L with R. Similar Hamiltonians, with appropriate changes of the quark fields 
and the CKM matrix elements describe also B®-B® and K°-K° transitions. In the SM 
only the Wilson coefficients of Q VLL are relevant and the details of their contributions are 
well known (see e.g. [37]). 



The basic formula for the B®-B® mass difference reads 



AM, 



r 2 m 2 

U F 1VI W M p2 n |T/eff|2| j?q\ 

6?r 2 M B q VB-r Bq BB q \Vt q I \*tt\ 



(6.6) 



where q = s or d. We recall also that [23] 
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with r = 0.985 [38] describing 0(a s ) QCD corrections to ^(xt) in the SM and t]b = 0.55. 
The factor f s introduced in eq. (6.3) is given by 1 + f s = Ff t /So(x t ). Next 

pa 



pa 



(6.8) 



^VbB Bs 

where B Bs is the non-perturbative parameter related to the matrix elements of <5 VLL . The 
coefficients P" include NLO QCD renormalization group factors [39, 40, 41] that sum up 
large logarithms between fi S and ^ = 0{m h ) scales. Explicit formulae for Pf in terms of 
these QCD factors and non-perturbative parameters Bi can be found in [41, 23]. 

In [23], except for B VLL , we have set the non-perturbative parameters Bi in the B- 
system to unity as the results of lattice calculations were not available at that time. 
Meanwhile all these parameters have been calculated in [42]. Translating these results 
into our operator basis by means of the formulae given in [41] we find 



b vll = ag7j b l,k = L75> b lk = L16j B ^ = 80j B suu = o 71 



->LR 



■>LR 



SLL 



?SLL 



(6.9) 



for /ib — 4.6 GeV. This gives 

P 1 Li? = -2.03, P 2 Li? = 2.56, P 1 5iL = -1.06, P 2 5LL = -2.05 (6.10) 
and using (6.8) with B Bs = 1.3 we find 



p£ R = -0.71, P 



LR 



0.90, P{ 



SLL 



-0.37, P 



SLL 



— 0.72. 



(6.11) 



We observe substantial suppression of APf LL relative to the vacuum insertion estimates 
used in [23, 27] that resulted in —0.53. The formula for e K in terms of Ff t as well as the 
expressions analogous to (6.7) for Fjr t can be found in ref. [23]. 
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In our numerical analysis of section 7 we have included complete expressions for all the 
one loop box diagrams contributing to Wilson coefficients Ci(fis) relevant for B° sd -B° sd 
mass difference and Ek as well as the complete expressions for the contribution of the 
double penguin diagrams shown in fig. 14. In the box diagrams involving charged Higgs 
and Goldstone bosons we have used vertices with corrections ej/^ and e^j^ calculated 
as in sec. 2. Similarly, in the diagrams shown in fig. 14 we have used full flavour changing 



vertices 



X RL 



ji 



X LR 



JI 



as given in eqs. (2.34), (2.35), (2.37) and (2.38). The necessary 
formulae are collected in the Appendix A. In the numerical calculation for the scale fis 
we have taken /is = M H +. 

Let us now concentrate on the most important of new contributions to the Wilson 



coefficients relevant for AM S . 



6.1.2 Double Penguin Diagrams 

As shown in [19, 20, 23] in supersymmetry with large tan (3 contribution of the double- 



iSLL 



penguin diagrams shown in fig. 14 to the Wilson coefficients of the operators Q% , Q\ 
and Qf RR can be substantial. Their Wilson coefficients Ci relevant for the B®-B° s mixing 
are then determined from the matching conditions: 



167T 2 
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tb*-\rts\ r<LR _ ■ 
cS V cS j U 2 — 1 
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X 



ibs 



RL 
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LR 
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(6.12) 



S°=h°,H°,A° 

with X^ L X[ R replaced by X| L X| L and X[ R X[ R for Cf LL and Cf RR , respectively. In 
these cases a combinatorial factor 1/2 has to be included on the r.h.s of (6.12). The 
contribution of the neutral Goldstone boson G° can be neglected as X^ R and X RL vanish 
in the SU(2) x U(l) symmetry limit. (See eqs. (3.23) and (3.24).) 

In what follows it will be convenient to introduce the "reduced" couplings 
defined by 



X s 
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RL 



bs 



ee Gf2'^V cS 

bs 



III* 



ts 



X LR 
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RL 



bs 



(6.13) 



with the same definition for 

In the approximation of sec. 3, using (3.26) and (3.27) we find for J 

WibfnfeY tan/5 



3 = b, I = 2 



X s 



'1 + e 3 tan/3)(l + e tan/5) 



(xf -xftan/3) 



and 



X 



LR 



bs 



m s m^€Y tan (3 



(xf -xf tan/3) 



(6.14) 



(6.15) 



(1 + e 3 tan (3) (1 + e tan (5) 

where we used the approximations 1/ cos/3 tan/3, 2\/2GF7nf valid for tan/5 3> 1. 
xf and xf are defined in (3.21). Following the numerical analysis performed in sec. 5.2 
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reliable approximation to the results of the complete calculation requires, however, using 
the formulae of section 3.6 with the epsilon parameters extracted as in eq. (5.1). Still, 
the simple formulae (5.2) used in (6.14), (6.15) describe well the qualitative behaviour of 
the double penguin contribution. 

From (6.12) and (6.13) we find then the double penguin contributions to the Wilson 
coefficients in question: 



and 



qSLL _ 



G F (16tt 
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lG F (167r 2 ) 2 v _L_ 



X RL 
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S =h°,H°,A<> 1V1 S° 

with Cf RR obtained from Cf LL by interchanging L and R 
Using formulae (6.14) and (6.15) we find 

mi>m 2 eY tan 2 j3 



bs r ^ 
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bs 



(6.16) 



(6.17) 
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m s m 2 ey tan 2 (3 



[sin(ct — (3), cos(a — j3), —i] 
sin(a; — (3), cos(a; — (3), i] 



(6.18) 



(6.19) 



'1 + e 3 tan/3)(l + e tan/3) 

for S = [H°, h°, A ]. Consequently the dominant, C(tan 4 /3), contributions to C^ R , Cf LL 
and C]" RR , that come solely from the double penguin diagrams read 

tan 4 /3 e y (167r 2 ) 2 



C 9 LR 



GFm b m s Wq 
V^^M 2 ? (l + e 3 tan/3) 2 (l + e tan/?) i 



(6.20) 
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(6.21) 
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s ^ySLL 
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(6.22) 



(Note that ey as given by eq. (5.2) is equal Xc t /167r 2 of ref. [23] and Gf/\/2 = 
7tq; em /2s 2 ^M^.) Here 



sin 2 (o; — pi) cos 2 (o; — P) 1 



M 2 



+ 



Ml, 



± 



M 2 M| M 2 



(6.23) 



where we have used the fact that in the MSSM for tan P >> 1 and M40 > M z one has 
sin 2 (a; — p) 1 and cos 2 (a — /3) 0. Because M40 ps M#o, it follows, that for large 
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tan f3 the H° and A contributions to Cf LL and Cf RR cancel each other [20] and the 
contribution of h° can be neglected. It is therefore the coefficient C L R that receives the 
largest contribution from double penguin diagrams [23] . Their contribution to the Wilson 
coefficients relevant for the B d -B d mixing (for ek) are given by the same formulae with 
m s replaced by (m s replaced by and replaced by m s ). Our result for C^ K 
agrees with the corrected version of eq. (42) of [27] where the presence of additional tan (3 
factors involving e 3 and e , not included in [23] has been pointed out. 

An important feature of the double penguin contribution to C^ R , (and also to Cf LL 
and Cf RR ) is its fixed negative sign that is the same as the sign of the contribution of the 
charged Higgs boson box diagrams at large tan (3 [23]. Because the strong correction factor 
P 2 LR in (6.7) is positive the double penguin contribution interferes destructively with the 
SM ones and leads to 1 + f s < 1. Another interesting feature is its strong dependence on 
the left-right mixing of the top squarks as ey oc A(. 

Using the general formulae given in refs. [23, 41] and collected in the preceding sub- 
section we find the contribution of the double penguin diagrams to AM,: 



(AM.) 



DP 



24tt 2 



B.Ft.\V*\ 2 P 2 LR CF(ns). 



(6.24) 



For us — 0(m t ) the factor P 2 LjR summarizing renormalization group effects between ^ < 
pL < /is an d including the relevant hadronic matrix elements 

Identifying C\ K in (6.20) with C~ r (ii s ) we find for large tan/5 
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3.0GeV 



m s (fi s ) 
0.06GeV 



2.50 



B s 



230 MeV 



Is 
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(167r 2 ) 2 e 2 Y 



(1 + e 3 tan/?) 2 (l +e tan/3) 5 



(6.25) 



6.1.3 Charged Higgs Box Diagram Contributions 

The formulae for the charged Higgs boson (H ± , G ± ) box diagram contributions to the 
Wilson coefficients of the effective Hamiltonians describing B Q sd -B G sd mixing have been 
presented in section 4.1 of [23] and have been confirmed in [27]. They do not, however 
take into account the corrections to the H + couplings discussed in sections 2.5 and 3.4. 
The corrections 6jf and 6jP can be easily incorporated in those expressions of eqs. (4.4) 
and (4.5) of [23] that contain a single D or a single _D 2 -function. To this end it suffices 
to make the following replacements 



Jik=i 

a L(R) 



pH 

r RL(LR) 



.11 



Mi 



„JIk=2 
L(R) 



pG 



JI 



/Kff (6-26) 

The two expressions in eqs. (4.4) and (4.5) of [23] which contain three -D 2 -functions have 
been obtained by using the unitarity of the CKM matrix. Therefore, in this case the 
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replacement (6.26) has to be supplemented by multiplication of the second and third D2- 
function in the square bracket in (4.4) of [23] by r and r 2 , respectively and by r 2 and r 4 
in (4.5) of [23] where r = (1 + €3tan/3)/(l + eitan/3) as follows from the formula (3.47). 
For completeness we give the relevant expressions in the Appendix A. 4. 

In this section we would like to asses the magnitude of the charged Higgs boson box 
diagram contribution to AM S and AM d with the corrections ejj^ included (as follows 
from the numerical analysis of sec. 5, the impact of the corrections e]j^ is negligible) 
and compare it with the magnitude of the contribution of the double penguin diagrams. 

As found in [23] the largest contribution to AM S of the charged Higgs box diagrams, 
arises from the coefficient C% R . This contribution is proportional to m s m(, tan 2 (5 and its 
sign is opposite to the SM contribution. The analogous contributions to AM d and e K are 
much smaller being proportional to m d m,b tan 2 (3 and m d m s tan 2 f3, respectively. 
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Figure 15: Comparison of the double penguin and charged Higgs boson box diagram 
contributions to the ratios AM S /(AM S ) SM (panel a) and AM d /(AM d ) SM (panel b) for 
Ma = 200 GeV and tan (3 = 50. Only points giving acceptable BR(B — > X s ^) are shown. 
Dashed lines show the ratios (AM S ) H+ /(AM S ) SM and (AM d ) H+ /(AM d ) SM in the absence 
of radiative corrections to the charged Higgs boson vertices. 

An order of magnitude estimate of the charged Higgs boson box diagram contribution 
to AM S for M H + PS m t can be obtained by taking into account only the contribution to 
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C^ R in (6-7) of the box diagram with two H ± : 

w „« m ( , 27) 

where P^ R ~ 2.56 and the factor (1 + e 3 tan/3) 2 in the denominator comes from the rule 
(3.47). One should remember, however, that the box diagram with the W ± H^ exchange 
also gives the contribution oc tan 2 (3 and with the same sign as the one with H ± H T 
exchange [23]. 

In figure 15 we compare for our scan over the MSSM parameter space the relative 
magnitude of the (complete) charged Higgs boson contributions to AM S and AM d with 
the contributions of the double penguin diagrams. We observe that in the case of AM S 
the contributions of the box-diagrams is much smaller than that of the double penguin 
diagrams for most of the points in the parameter space. The latter contribution is always 
negative and in a large portion of the parameter space exceeds also the SM one. The upper 
(lower) branch of points seen in panel a of figure 15a corresponds to \x > (pi < 0). In 
fact, for points corresponding to unnaturally large A t combined with negative values of fi 
leading to enhancement of the effective couplings, the ratio (AM S ) DP /(AM S ) SM can reach 
values even as big as —10. The same points lead also to (AM d ) DP /(AM d ) SM < —0.2. They 
are excluded as they lead to unacceptably big BR(B® — > /i + /U~) - see sec. 6.3. In contrast, 
the double penguin contribution to AM d can be also positive. This is because in this case 
the suppression by rn^frid of the Wilson coefficient C\ R is stronger than the suppression of 
Cf hL by the factor T~ in eq. (6.21). Since P^ LL in (6.7) is negative this leads to positive 
contribution to (AM d ) DP /(AM d ) SM . As is seen from fig. 15 for some special values of the 
MSSM parameters (AMd) DP /(AMd) SM can reach ~ ±0.2 but is much smaller in most 
of the parameter space. The contribution of the box diagrams involving H ± to AM^ for 
large tan j3 is always negligible compared to the SM and double penguin contributions. 

6.1.4 Chargino Box Diagram Contributions 

In the scenario considered in [23] and in the present paper supersymmetric particles are 

heavier than the Higgs bosons and the chargino box contribution (AM S ) X± is small as it 

is suppressed by 1/M| USY . We have checked that in our scan, with all sparticles heavier 

than 500 GeV, the ratio (AM S ) X± /(AM S ) SM can reach only 0.02 for acceptable values of 

the B — > Xg'-f rate. However, we find some differences with the formulae derived in [27] 

and we want to devote this subsection to clarify this issue. 8 

8 P.H.Ch. would like to thank Gino Isidori for the discussion which allowed to identify the source of 
the discrepancy. 
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The contributions of box diagrams with chargino exchanges to AM Sjd at large tan (3 
have been calculated in [23, 27] but only in [27] they have been presented in details. These 
authors have found that the most important contribution of these diagrams is present in 
the coefficient Cf LL that behaves as m^m 2 tan 2 (5 and contributes to AM S](i with the same 
sign as the SM contribution. As stressed in [27], if the sparticles are light this contribution 
could be relevant also for AM d due to the absence of the suppression factor m d . These 
statements have to be revised, however. 

First of all, let us notice that, as follows from the exact formulae for the chargino 
contributions collected in the Appendix A. 4, the Wilson coefficient Cf LL in (6.4) has 
exactly the same dependence on the chargino couplings as does Cf LL . In fact, the term in 
the chargino box amplitude which has that combination of the chargino couplings gives rise 
to a single scalar left-left operator which, however, has colour indices contracted differently 
than the operators in the basis (6.5) used also in [27]. Making the Fierz transformation of 
this operator gives rise to the Of LL and (9f LL operators with the Wilson coefficients given 
in the Appendix A. 4. The authors of ref. [27] incorrectly identified the scalar left-left 
operator with wrong contraction of the colour indices with Of LL . This has the following 
consequences. Firstly, our Cf LL has an additional factor —1/2 which results from the 
Fierz transformation. This reverses the sign of the contribution of Cf LL and decreases it 
by one half as compared to ref. [27]. Secondly, C| LL = — (l/4)Cf LL also contributes to 
the mass difference and further reduces the effects of Cf LL . More precisely, using (6.7) 
and (6.8) the tan 2 /? contribution of charginos to (6.6) enters in the combination 

where we have used the fact that P 2 SLL « 2P X SLL (see eq. (6.10)). This reduces the tan 2 /3 
contribution of charginos by roughly another 1/2. 

Thus, our dominant at large tan/? part of the chargino contribution to AM Sjd is ef- 
fectively smaller by a factor of approximately 4 compared to the one of ref. [27] and has 
opposite sign compared to the SM contribution. Neglecting the gauge coupling constants 
as well as the Yukawa couplings other than top and bottom one gets for the dominant 
chargino box diagram contribution 

r SLL d, \- o m l m t A2 t tan 2 /3 2 2i 7 2i 7 2j 7 ij 
° l " ~ 2 (1 + 6 3 tan^) 2 £ Z ~ ^ Z ' 

x [A)(l, 1) - 2D (1, 2) + D (2, 2)] (6.28) 

where t 1 is the heavier stop and D (l, k) stands for D (m Cl , mc d , M^ k , Mf ). The epsilon 
factor in the denominator appears here in agreement with the rule (3.48). With all these 
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Figure 16: Comparison of the double penguin and chargino box diagram contributions 
to the ratios AM S /(AM S ) SM (panel a) and AM d /(AM d ) SM (panel b) for M A = 200 GeV 
and tan (3 = 50. Only points giving acceptable BR(B — > X s 7) are shown. In this scan we 
allowed for charginos and stop as light as 100 GeV. 

corrections we find that even for tan (3 ~ 50 the contribution of charginos to AM s ^ d is in 
general dominated by the C VLL Wilson coefficient which has the same sign as the SM. 
Only for very big values of A t can the chargino contribution to the Wilson coefficients of 
the scalar operators compete with the contribution of the standard VLL one. 

In figure 16 we show the relative magnitude of the chargino box and double penguin 
contributions to AM S and AM d for a scan which included charginos and stop as light as 
100 GeV. The biggest contribution of charginos allowed by B — * X s ^y constraint is ~ 0.2 
of the SM one. If the B -> X,7 constraint is relaxed, the ratios (AM S ) X± /(AM S ) SM and 
(AM d ) x± / (AM d ) SM can reach ~ 0.3. The anti-correlation of the chargino and double 
penguin contribution is clearly seen in figure 16 (it is even better visible if the points with 
unacceptable BR(B — > X s ^) are retained). The reason for it is simple: the biggest double 
penguin contribution arises for small values of (1 + e3tan/3) when also the contribution 
the scalar left-left operators generated by boxes is enhanced and cancels the contribution 
of the O vll generated by the chargino box diagrams. The biggest total chargino box 
contribution arises when (1 + e 3 tan/9) is as big as possible to suppress their negative 
contribution through the scalar left-left operators. For ji < this requires A t > which is 
incompatible with BR(B — > X s '~f) constraint. For // > the biggest values of (l + e 3 tan/3) 
are obtained for A t < which again leads to unacceptable B — > X s ^ rate. (Cf. the signs 
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of ey and eo in eqs. (5.2).) 



6.2 The B% 

s,d 



H + H Decay 



The effective Hamiltonians describing these decays has the form 
2G F 
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eff 



2irsm 2 9 w 
where q = s or d and 



V t f*V t f[c A A + c' a O' a + c s Os + c' s O' s + c p O p + c' p O' p ], (6.29) 



O s = m b (b^q L )(ll), 0' s = m q (hq R )(ll), 
Op = m b (b R q L )(l'y 5 l), O'p = m q (b L q R )(l^ 5 l). 



(6.30) 
(6.31) 
(6.32) 



For large tan/5 the contributions from neutral Higgs scalars dominate. Using the 
Lagrangian (2.33) for the q = s case and neglecting m s with respect to rrib we find 
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with the reduced coupling 
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defined in (6.13) and 
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Comparing (6.33) with (6.29) and using (6.14) we find 
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(6.36) 



The contributions of the operators 0' s and O'p are strongly suppressed by m s with 
respect to Os and Op. Consequently following [43] we can write to an excellent approxi- 
mation 
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where cs and cp are the dimensionless Wilson coefficients 



c s = M Ba c s , c P = M Bs c P . 



(6.38) 



The coefficients ca and c' A receive contributions from Z°-penguin diagrams and box dia- 
grams and are weighted by the muon mass. For large tan (3 they can be safely neglected 
with respect to cs and dp that grow like tan 3 /3 and are not chirally suppressed by the 
muon mass. 

In the large tan (3 limit the contribution of h° to cs can be neglected and setting 
M%> we find from (6.35) and (6.36) 



M 2 H0 
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with cp given in (6.36). Consequently we find the branching ratio 
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This result agrees with [27]. Moreover one has 
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(6.41) 



that is, the ratio of the branching fractions can depend on the SUSY parameters only 
weakly through \V^f /V^f\ which should be consistently determined from the unitarity 
triangle analysis. 

The presence of additional tan (3 dependence in the denominators of eqs. (6.25) and (6.40), 
not included in [23] and [20, 6, 22, 21], has been pointed out in [27]. While we confirm 
these additional factors, we would like to emphasize that depending on the sign of the 
super symmetric parameter \x they can suppress AM® P and BR(B° S — > relative to 

the estimates in the papers in question, as stressed in [27], but can also provide additional 
enhancements. 



6.3 Correlation between B° o , 

s,d 



and AM DP 



We are now in the position to give an explicit formula for a correlation between the neutral 
Higgs contributions to B Q S — > li+ it and AM® P that we have pointed out in [23, 29]. 
Indeed from (6.25) and (6.40) we find 
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where 
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(6.43) 



This relation depends sensitively on M A o and tan/5 but it does not depend on e and 
e 3 . From (6.41) a similar correlation between BR(B® — > yU + /i~) and AMf p follows. If 
the flavour dependence in the epsilon parameters is taken into account as in section 3.6 
and formulae (3.55) and (3.56) are used, the r.h.s of (6.42) is multiplied by the factor 
e (32)^ e (23) r32 ^ -y^^j^^ however, does not depart significantly from unity. In our numerical 
analysis of section 7 we use the formulae of the full approach anyway. 

In order to understand these results better, let us now assume that AM S has been 
measured and that appropriate supersymmetric parameters can be found for which the 
MSSM considered here agrees with (AM s ) cxp . If < (l + / a ) < 1 this implies (AM s ) exp < 
(AM S ) SM . Then combining (6.3) and (6.42) we find [29] 
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(6.44) 



with "=)=" corresponding to < (1 
analogous expression for BR(B® — 



f- f s ) < 1 and (1 + f s ) < 0, respectively. Using (6.41) 
can be found. In writing (6.44) we have taken 
into account that (AM S ) DP is always negative and that for large tan f3 (AM S ) H± is negative 
and (AM,) X± mostly positive. Formula (6.44) is valid provided the expression in square 
brackets is positive and larger than 10~ 3 . Otherwise, other contributions, in particular 
those coming from Z°-penguins have to be taken into account. In our numerical analysis 
we take them into account anyway. 

Formula (6.44) demonstrates very clearly that if (AM s ) cxp will turn out to be close or 
larger than the SM value, the order of magnitude enhancements of BR(B Q s d — > /i + /i~) in 
the scenario of the MSSM considered here with < (1 + f s ) < 1 will be excluded. On 
the other hand large enhancements of BR(B° s d — > /x + /x~) are in principle still possible if 
the double-penguin contribution is so large that (1 + f s ) < and the "+" sign in (6.44) 
applies. For tan/3 < 50 obtaining (1 + f s ) < and the right magnitude of AM, requires 
fi < so that the couplings (3.26) are enhanced by the e-factors in the denominator, 
/i < is excluded in particular scenarios like minimal SUGRA, in which the sign of A t is 
fixed and /i < does not allow for satisfying the B — > X s ^y constraint [26], but cannot be 
excluded in general. 

We will analyze (6.44) numerically in detail in section 7.3. 
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6.4 Unitarity Triangle 

The rate of the B d — > decay depends on the element V t e f of the low energy CKM 

matrix. The value of V t e f cannot be extracted by using exclusively tree level dominated 
processes and requires an analysis of the combination of data on \V*f? /V^ b |, sin2/? ut from 
the CP asymmetry a^^ s , Ek, AM d and AM S (the so-called UT analysis). While in 
the scenario of the MSSM we are considering the extraction of the first two quantities 
from the data is not affected by the new particles, Ek, AM d and AM S receive additional 
contributions which can modify the extracted value of V t °f . The unitarity triangle in the 
MSSM with large tan/? has been discussed in refs. [23, 7]. Here we recall this discussion 
briefly and outline the procedure we use in sec. 7.3 in order to consistently determine V t e d s 
for a given set of the MSSM parameters. 

Let us first recall that the contribution of the new particles enter ek, AM^ and AM S 
through the factors 



Ftt = S (x t )[l + f £ ], F* = S (x t )[l + f d ], F t s t = S (x t )[l + f s ] (6.45) 



respectively where So(x t ) is the universal contribution of the SM particles. F t s t is explicitly 
given in eq. (6.7). The hierarchy of the new contributions to the parameters f\ is then as 



• The contributions to f £ at large tan (3 generated by the diagrams similar to the 
ones shown in figure 14 are fully negligible being strongly suppressed by m s m d , m 2 s 
and higher order terms in m d)S . Also the charged Higgs and chargino box diagram 
contributions are known to be negligible for tan/3 ^> 1. Consequently f e can be set 
to zero and the constraint on the UT stemming from Ek is the same as in the SM. 

• The double penguin and charged Higgs boson contributions to AM d are proportional 
to mdfrib tan 4 (3 and m d fnb tan 2 f3, respectively [23, 27]. Detailed calculation shows 
that \f d \ is at least one order of magnitude smaller than |/ s |. The latter, in the 
most likely scenario with < 1 + f s < 1, is in turn bounded by (AM s ) exp to be 
smaller than 0.5 [23]. Consequently \f d \ ^ 0.05 and can be safely neglected in view 
of the uncertainties in F Bi \J B Bd . In the unlikely scenario in which 1 + f s < and 
\f s \ > 1-5, \fd\ could reach 0.2 and having negative sign would suppress AM d relative 



to the SM prediction. For fixed F Bd \J B Bd and (AM d ) cxp , this would imply a larger 



• The new contributions to f s can be substantial [23, 29]. They will be investigated 
numerically in detail in section 7 using the approach developed in sec. 2 which takes 
into account the resumation of leading higher order terms pointed out in [27]. In our 



follows: 
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scenario f s is always negative and suppresses AM S with respect to the expectations 
based on the SM. 

In order to determine the impact of f s ^ on the determination of \V t e f\ let us recall 
that 

AM d ^F 2 B B Bd \V t f\ 2 S (x t ) (6.46) 
AM S oc F 2 B B Bs \V t f\ 2 S (x t )\l + f s \ (6.47) 
and 




(6.48) 



Exact expressions can be found in [23, 37]. 

If the value of F Bd Bs d was known precisely \ V t c f | could be determined directly from eq. 
(6.46) and the value of f s would not have any impact on it. In our scenario \V t e d -\ would 
then assume the same value as in the SM. Since |V^ ff | is almost fixed by the unitarity of 
CKM matrix the eq. (6.47) would then only provide a constraint on the MSSM parameters 
space, the more stringent, the smaller was the uncertainty in F Bs Bb s [23]. In other words, 
a precise value of F Bs Bb s combined with a precise measurement of AM S could directly 
probe by means of the formula (6.47) the value of |1 + / s |. 

At present, however, F Bd Bs d (and also F Bs Bb s ) is poorly known and a more precise 
value of \V td s \ is obtained by using the relation (6.48) in which £, calculated by QCD sum 
rules or lattice methods that are insensitive to new short distance contributions, is known 
with much better accuracy than are F Bd Bs d and F Bs Bs a separately. This means that 
F Bd B Bd and F Bs B Ba are positively correlated and cannot simultaneously assume values 
from the opposite extremes of their respective allowed ranges. Therefore, as seen in (6.48) 
< 1 + f s < 1 leads to a smaller value of \V t e d | and smaller angle 7 than does f s = 0. For 
example, if the value of AM, is close to its present experimental bound, |1 + f s \ 0.6 
gives \ V t f\ « 7x 1CT 3 . This can be also viewed as follows: for a given experimental value 
of AM S , the factor F^ g B Bs must be larger if < |1 + / s | < 1 than if f s = 0. With £ fixed, 
this means that also F Bd B Bd must be larger implying a smaller value of \V t e f\ and of the 
angle 7. Similarly, 1 + f s < leads to bigger \V t e d s \ and bigger 7 provided |1 + / s | > 1 For 
example if 1 + f s m -1.3 one has \V t f\ w 1 x KT 2 [7]. 

In our global scans in sec. 7 we proceed as in [29, 7]: For a given set of the MSSM 
parameters we scan over the Wolfenstein parameters A, A, p and fj which parametrize 
V eS (see e.g. [37]) as well as over the non-perturbative parameters F B< B Bd and F\B Bs in 
their respective ranges specified in ref. [44] and compute the quantities of interest only 
for those A, A, p and f\ for which e K) AM d) a^K s = sin-2(3ut, iKf/Klfl an d £ assume 
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acceptable values. This allows to take into account in BR(B® — > fi + fi ) also the small 
effects related to the dependence of |V^ ff | on the supersymmetric parameters. 

6.5 B -> X s j 

Finally we recall certain aspects of the BR(B — > X s ^) calculation referring frequently 
to the formulae present in the literature. We would like to present for completeness 
simple recipes allowing to use the corrections e HL etc. introduced in sec. 2 in the existing 
formulae for the relevant Wilson coefficients. We also show how to include the recently 
found correction depending on the neutral Higgs boson masses [28] and on the couplings 
[Xlr] JI in these formulae. 

As far as the SM part of the computation is concerned, we closely follow the approach 
of ref. [45]. The two loop contribution of H + is taken from ref. [46] (see also [47]). As in 
ref. [25], for our scenario with heavy sparticles we include only one loop contribution of 
the supersymmetric particles as given in refs. [48, 33]. 

6.5.1 Charged Goldstone and H + Contributions 

In order to see how to include our corrections in the calculation of the B — > X s ^ rate let us 
recall that at large tan (3 the one-loop contributions of G + (H + ) to the Wilson coefficients 
are proportional to the product of the tLG + {H + )bn vertex and the Hermitian conju- 
gate of the tnG + (H + )sL vertex. Neglecting the contributions originating in the chiral flip 
on the external b quark line and not including yet the large tan/3 enhanced corrections 
one has 

CJ<? = -Fi%(x t ), = F^(y t ) . (6.49) 

Here x t = rn^/M^, y t = ml/Mf I+ . Combining the G + and W + contributions results in 
replacing —F^(x t ) by F^(x t ). The explicit expressions for F^(x t ) and F^(x t ) can 
be found in [25]. 

Let us consider the G + contribution first. As we have demonstrated in sec. 3, in the 
577(2) x U(l) symmetry limit ejf = e^f = and the couplings of G + have the same 
form (3.38) as in the SM. This means that in the SU(2) x U(l) symmetry limit this 
contribution is correctly accounted for in the SM contribution to C$L. Small departures 
of ej/^ and ejf from zero which arise beyond this approximation can be also included but 
our numerical studies of sec. 5 show that they are too small to have any impact on the 
B — > X s 7 rate. Nevertheless, we can include them by noting that the full contribution 
G + to the Wilson coefficients is 

Ci? 8 V = -(1 + e t G /)(l + 4 R )F%hxt). (6.50) 
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As for ef s L = ef b R = this contribution is taken already into account in the SM contribu- 
tion, one has only to modify by 

* G+ C% = ~(eg R + eg L + eg R ef s L )F^( Xt ) « -(eg* + egV$>(*«) (6.51) 

where we have neglected the higher order term. The prescription (6.51) replaces in our 
approach the recipe given in eq. (18) of ref. [25]. It formally agrees with the latter if 
one drops ef s L not considered in [25] and identifies Vt s in the inG + SL vertex used in [25] 
with V£f (see our discussion in section 4) and recalls that V# = V tb . Using then the 
substitution rule (4.4) in the opposite direction one finds 

- e GR _^ 1 _ l + 4(t)tan/3 = [e b -e' b (t)}t&n(3 

tb ' 1 + e 6 tan/3 1 + e fe tan/3 1 ' 1 

as in eq. (18) of [25]. We have stressed however, that one has to be careful in evaluating the 
corrections ef b R and ef s L , because uncontrolled approximation can result in overestimating 
their impact on the B — > X s 7 rate. 

In the case of H + contribution, the formula analogous to (6.50) reads 

<' = (l-fl(l-4VaW. (6-53) 
In the approximation of section 3, using eqs. (3.43) and (3.47), we have: 



+ eo tan (3 



(6.54) 



in agreement with the factor multiplying F^(y t ) in the formula (87) of [28]. However 
as shown in section 5, the approximation (3.44) is not numerically very accurate and 
better approximation is obtained by using for ef s L eq. (5.8) with J — 3, I — 2 and 
e tl R = e 3 tan/3/(l + tan/3) obtained from (3.59) with e 3 calculated by means of (5.1). 
Thus, the effect of corrections to the charged Higgs boson vertices results in the following 
modification of the Wilson coefficients c!$L : 

^cJg) = K L - 4* + 4 L 4 R ) F u8)(yt)- (e.ss) 

Sticking to the naive approximation of sec. 3 with the Ajj term in (3.43) neglected, i.e. 
inserting in (6.55) 

d L = e' t (s) tan/3, e» R = - - (6.56) 

•ef. [25] if, in agreemer 
identifies V ts with V t *f in the t R H + s L vertex of [25]. 



■is ~tw— -to l + ebtan (3 

one recovers the rule (19) of ref. [25] if, in agreement with our discussion of sec. 4, one 
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6.5.2 Chargino Contributions 

The inclusion of the tan (3 enhanced corrections in the chargino contributions to Wilson 
coefficients of higher dimensional operators involving down quark fields of arbitrary flavour 
is achieved by means of the rules (2.53) and (2.54) or (3.48). For the couplings $c 1 l sl 
(2.53) amounts simply to the replacements 

mt - m u V ts - V t f. (6.57) 

For the couplings fic l R bn the rule (2.54) or (3.48) applies implying 

where following the refinements of sec. 5 e 3 = + y\ 6y 3 ^ • This is accidentally equivalent 
to expressing in the tree level formulae simply V tb and m b through V t c b s and mj, respec- 
tively. The formula (30) of ref. [25] with e b = e ds = eg 3 -* + y t 2 ey ^ is consistent with this rule 
rule provided V ts and V t b in this formula are identified with the effective CKM elements. 

On the other hand for the Wilson coefficient of the m s /m b suppressed operator m s s R a^ u b L 
the situation is more involved as our rules imply 

m -> rfk, V tb -> V t f (6.59) 
for the vertex involving the quark and 

for the vertex involving the quark. If the flavour dependent effects are important 
1/(1 + e 3 tan/5) in (6.60) should be replaced by the square bracket of the formula (3.61) 
with 1 = 2. Note that expressing in the tree level formulae simply V ts and m s through 
and m s , respectively would give wrong tan/3 dependence. 

6.5.3 Neutral Higgs boson contribution 

For completeness we include the neutral Higgs contribution to the Wilson coefficients 
found in ref. [28]. It is the natural consequence of the flavour changing neutral Higgs 
coupling [X L ^] 23 which, through the diagram shown in figure 17, leads to the contribution 
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Figure 17: Additional diagram contributing to b — > 57 or b — ■> transition. The dot and 



the cross denote the flavour changing coupling [XLn] sb and the helicity flip, respectively. 



where 



33 



and 



23 



are given in eqs. (3.22) and (3.27), respectively. 
The last comment concerns the dependence of the B — > X s 7 amplitude on the CKM 
matrix elements. In our scans presented in sec. 7 for a given set of the MSSM parameters 
we use the value of \V^ S determined consistently from the UT analysis as described 

in the preceding subsection. We do not need therefore to implement the recipe of ref. 
for correcting the B — > amplitude for new physics effects in V^f . 



7 Numerical Analysis 

In this section we present numerical analysis of the dependence of AMd, s and B® d — > 
on the parameters of the MSSM. We will also present the global analysis of these quantities 
taking into account available experimental constraints, in particular the one from the 
measured rate of the B — > X s ^ decay. We present the results based on our complete 
approach of section 2 which includes automatically the £77(2) x U(l) breaking corrections 
as well as the dependence of the flavour changing couplings on the electroweak gauge 
couplings. On some plots we compare these results with the one obtained by using the 
approximation of sec. 3 based on SU(2) x £7(1) limit and dominance of a s and the top 
and bottom Yukawa couplings. The latter describe qualitatively the main features of the 
MSSM effects but are not very accurate. 

7.1 The Size of f s and f d 

The parameters f s and fd introduced in eq. (6.3) and directly related to the ratio 
AM M /(AM M ) SM : 

AM M /(AM M ) SM = |1+/ M | 

receive contributions from double penguins, charged Higgs boson box diagrams and chargino 
box diagrams. As we have already said in sec. 6.1.4, for sparticles heavier than 500 GeV 
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the latter contributions can be safely neglected. We have also shown in sec. 6.1.3 that for 
large tan (3 the contribution of the H ± box diagrams to f s is negative and can give up to 
S H f s ~ —0.1 for /i > and up to S H f s ~ —0.2 for jj, < 0. Except for very small values of 
the At parameter, for which the double penguin diagram contribution is small, the charged 
Higgs boson box diagram contributions are much smaller than the latter one and do not 
play any role in the correlation between AM S and BR(B° s d — > The contribution 

of the charged Higgs boson box diagrams to f d is negligible —0.02 < 5 H f d < whereas 
the contribution of the double penguin diagrams to fd can be of either sign and for very 
special values of the parameters can reach \5 DF fd\ ~ 0.2. 




©t(deg) t (deg) 



Figure 18: (1 + / s ) in the MSSM for tan (5 = 50, M H+ = 200 GeV, m- g = 3M 2 , M 1r = 800 
GeV and the lighter chargino mass 600 GeV as a function of the stop mixing angle 6 t . In 
panel a) /i < and M h = 600 GeV, M h = 750 GeV. In panel b) fi > and M h = 500 
GeV, Mi = 850 GeV. Solid lines correspond to the complete calculation. Dotted ones 
- to the approximation based on the formulae of sec. 6.1.2 and the dashed lines to the 
calculation without the resummation of the tan (3 enhanced terms. Consecutive solid and 
dotted lines correspond to M 2 /|/i| = 3/4, 1 and 5/4 counting from left to right (in both 
panels). The results without the resummation do not depend in a visible way on M 2 /|/x|. 

In figs. 18 and 19 we show the factor (1 + / s ) as a function of the parameters to which 
its value is most sensitive, that is M H +, tan/5 and the mixing angle of the top squarks 6 t . 
The sign of the latter is always such that fiA t is positive, allowing for the cancellation of 
the charged Higgs boson and chargino contributions to the amplitude of the B — > X s ^ 
decay. 

Strong dependence of the double penguin contribution to f s on the left-right mixing 



63 



of the top squarks is clearly visible in fig. 18 where we show 1 + f s as a function of the 
stop mixing angle 9 t . As expected from our considerations of sec. 5, for // < (// > 0) 
the resummation of the tan/3 enhanced terms leads to dramatic increase (decrease [27]) 
of \f s \ due to the appearance of the factors (l + e 3 tan/3)(l + eotan/3) in the denominators 
in eqs. (6.20), (6.21) as compared to the formulae (4.11) of ref. [23]. Curves for different 
values of the ratio M 2 /\/jl\ in fig. 18 illustrate the impact of the heavier chargino. Indeed, 
for shown values of M 2 /\/i\ the lighter chargino is always almost higgsino-like and to an 
accuracy of few GeV, mc x ~ \fJ>\- (Without the resummation, the lines corresponding to 
different values of M 2 /\n\ are indistinguishable on the plot.) 

Dotted lines in fig. 18 have been obtained by using the formulae (6.20)-(6.22) and 
correspond to the approximation of sec. 3.3. We observe that for /i > (which is favoured 
by models like minimal SUGRA) when the resummation suppresses the effects of the 
double penguin diagrams as compared to the calculation of [23] , the approximation of sec. 
3.3 works well for \0 t \ < 20° but in general gives too strong a suppression. However, for \i < 
0, the approximation of sec. 3.3 grossly overestimates the impact of the double penguin 
diagrams on AM, already for 9 t < 5° and should not be used for realistic applications. 




200 400 600 800 1000 10 20 30 40 50 60 



M„ + (GeV) tan/3 

Figure 19: (l + / s ) in the MSSM for the lighter chargino mass 750 GeV, |r| = M 2 /\fi\ = 1, 
rrig = 3M 2 and M~ bR = 800 GeV as a function of M#+ for tan/? = 50 (panel a) and as 
a function of tan/3 for Mu+ = 200 GeV (panel b). Solid and dashed (dotted and dot- 
dashed) lines correspond to M h = 500 GeV, M h = 850 GeV (600 and 750 GeV). Solid 
and dotted lines, correspond to \x < whereas the dashed and dot-dashed lines to /i > 0. 
The stop mixing angle 8 t = +(— )10° for /j, < (// > 0). 

For the same value of the mixing angle 6 t , larger effects are obtained for bigger stop 
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mass splitting because in this case the parameter \A t \ has to be larger. This is clearly 
seen in fig. 19. It should be also stressed that the double penguin contribution does not 
vanish when the mass scale of the sparticles is increased (i.e. when all mass parameters 
are scaled uniformly). Thus, large effects decreasing 1 + f s below unity can be present 
in the MSSM also for the heavy sparticles provided the mass scale of the MSSM Higgs 
sector remains low and tan/3 is large as illustrated in fig. 19. Positive contribution to 
1 + f s seen in fig. 19b for tan/3 < 2.5 and M H + = 200 GeV is due to the ordinary charged 
Higgs boson box diagrams which contribute to the universal part of f s , fd and f e through 
the Wilson coefficient of the standard QY LL operator. For lighter H + and light charginos 
1 + f s can reach values ~ 2 [4] . Such high values of 1 + f s will be definitely excluded if 
the measured AM S is below 40/ps. 

The lines corresponding to M~ tx = 500 GeV, M^ 2 = 850 GeV in fig. 19 can be directly 
compared with their counterparts in fig. 13 of ref. [23]. As can be read off, for /i = —750 
GeV (/i = 750 GeV) the resummation of the tan (3 enhanced terms increases (decreases) 
f s roughly by a factor of 1.5 (2.4). 
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Figure 20: The ratio BR(B° S -> fi + BR(B° S -> /i+/i") SM in the MSSM for tan/3 = 50, 
Mh+ = 200, rrig = 3M2, M~ b = 800 GeV and the lighter chargino mass 600 GeV as a 
function of the stop mixing angle 9 t . In panel a) n < and Mi = 600 GeV, Mi = 750 
GeV. In panel b) fi > and Mf = 500 GeV, M^ 2 = 850 GeV. Solid lines correspond to 
the complete calculation. Dotted ones - to the approximation of based on the formulae 
of sec. 6.2 and the dashed lines to the calculation without the resummation of the tan/3 
enhanced terms. Consecutive lines of each type correspond to M 2 /|//| = 3/4, 1 and 5/4 
counting from left to right (in both panels). 
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7.2 Predictions for B® , — > /Lt+ \i 

In [6] and [22] the branching ratios _Bj? d — > have been calculated diagrammatically 

with one loop accuracy. The resummation of large tan (3 effects have not been done in 
these papers. The latter effects have been included in the SU{2) x U(l) limit first in 
[20] and their importance has been subsequently emphasized in [27]. Here we present 
for the first time results for B° sd — > fi + /i~ in the approach that combines the full one 
loop diagrammatic calculation of refs. [6] and [22] with the resummation of large tan f3 
enhanced terms. This includes both, SU (2) x U(l) breaking corrections to the calculations 
of refs. [20, 27] and the effects of the electroweak gauge couplings. We will asses the 
importance of these improvements with respect to the approaches of refs. [20, 27]. 

In figure 20 we show the ratio of the branching ratios BR(B Q S — > predicted in 

the SM and in the MSSM with tan (3 = 50 and M H + = 200 GeV as a function of the 
mixing angle of the top squark. Strong dependence on the latter parameter is evident. 
As in figure 18, curves of a given type corresponding to smaller values of M 2 /\fi\ illustrate 
the effects of lowering the mass of the heavier chargino (m^ ~ |/z|, mc 2 ~ M 2 fc> r all 
these curves) and of the increasing deviation of the lighter chargino from pure higgsino. 
Dotted lines in figure 20 have been obtained by using the formulae (6.35) and (6.36) 
of sec. 6.2 which correspond to the approach of ref. [27] that is, to the one of sec. 3 
with e and ey computed from the formulae (5.2) and (5.5), respectively. Dashed lines 
show the prediction of the MSSM without the resummation of the tan (3 enhanced terms. 
SU(2) x U(l) symmetry limit (sees. 3.2-3.4) of ref. [27]. 

As it could be expected, for /i > (/i < 0) the resummation decreases (increases) 
the branching ratio predicted in the MSSM. Similarly as in the case of f s , for /x > the 
approximation of sec. 3 gives the branching ratio smaller than the complete calculation, 
while for /i < it overestimates this rate significantly. For example, for // > 0, M 2 /|//| = 1 
and 9 t = —15° and other parametrs as in fig. 20, BR{B® — > /U + yU~) obtained in the 
approximation of sec. 6.2 is smaller than the one obtained using the complete approach 
by a factor of 1.5, while it is bigger by a factor of 2 if no resummation of tan/3 enhanced 
terms is performed. For /i < and 9 t = 15°, instead, the formulae of sec. 6.2 give 
BR(B® — > /U + /i~) approximately 3 times bigger than the complete calculation, while 
neglecting the resummation results in the branching ratio that is smaller by a factor of 

1.3 than obtained in the full approach. 

In figures 21a and 21b we show the dependence of the ratio of the branching ratios 
BR(Bg — > /i + /U~) predicted in the SM and in the MSSM as a functions of the charged 
Higgs boson mass and tan /3, respectively for two different choices of the top squark masses 
and both signs of \i. Here we show only the results of our complete calculation based on 
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Figure 21: The ratio BR(B° S -> fi+fi-)/BR(B° -> / u + / u~) SM in the MSSM for the lighter 
chargino mass 750 GeV, |r| = M 2 /\fi\ = 1, = 3M 2 and = 800 GeV as a function 
of M H + for tan/? = 50 (panel a) and as a function of tan/3 for M H + = 200 GeV (panel b). 
Solid and dashed (dotted and dot-dashed) lines correspond to = 500 GeV, Mf 2 = 850 
GeV (600 and 750 GeV). Solid and dotted lines, correspond to /j, < whereas the dashed 
and dot-dashed lines to /i > 0. The stop mixing angle 9 t = +(— )10° for /i < (/i > 0). 

the approach of sec. 2. 

The values of the ratio BR(B° S -> [i + / BR(B° S -> //+//~) SM shown in figures 20 
and 21 are also representative for the ratio BR(B® — > fi + BR(B® — > /i + /i~) SM if one 
neglects the small variation of |(V^f ) MSSM /(V^f ) SM | 2 with the supersymmetric parameters 
which we have discussed in sec. 6.4. 

Finally, in figure 22 we show the scatter plot of the ratio of the B® — * rates 
computed by using the formulae (6.35) and (6.36) of sec. 6.2 and computed using our 
complete approach of sec. 2 for the charged Higgs boson mass equal 300 and 400 GeV. 
Large, up to 50%, deviations of the approximate method with respect to the full calcu- 
lation are typical for most of the points. We observe that the B — > X s ^ constraint allow 
the approximate predictions of the branching ratio to be both, smaller and bigger than 
the one based on the complete calculation. Approximate predictions overestimating the 
exact ones correspond to very big negative contributions to AM S (1 + f s < 0) but, as 
found in [29], are not excluded by the lower experimental limit AM, > 15/ps provided 
1 + fs ~ —0.5. However, for M H+ = 200 GeV such points give BR(B® — > /i + /i~) above 
the CDF bound (6.2). 
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Figure 22: The ratio BR(B G S 



— \ approx 



>BR(B° S 



\ exact 



in the MSSM for 



tan/3 = 50 and Mh+ = 300 GeV (panel a) and 400 GeV (panel b). Points corresponding 
to unacceptable B X s ^ rate have been rejected. 



7.3 Correlation between AAf s and d 



As seen from eq. (6.44), important in the study of the correlation in question is the 
ratio (AM s ) exp /(AM s ) SM [29]. In order to find (AM s ) cxp /(AM s ) SM one has to deal with 
the non-perturbative uncertainties contained in the evaluation of (AM S ) SM . The allowed 
range for (AM s ) cxp /(AM s ) SM can be obtained by varying all relevant SM parameters like 
m t , V t f and Fb s \[Bb s - A conservative scanning of these parameters performed in [23] 
resulted in 



(AM S 



,exp 



15/ps 



< ^ M ^Z < b 



(AM, 



cxp 



15/ps 



(7.62) 



(AM S ) SM - 

with a = 0.52 and b = 1.29. It is however clear that the numerical values of the parameters 
a and b depend on the error analysis and the difference b — a should become smaller as 
the uncertainties in the parameters m t , and in particular in F Bs \J Bb 3 are reduced 
with time. For example, the very recent analysis which uses the Bayesian approach gives 
a = 0.71 and b = 1.0 [44] that correspond to the 95% probability range 15.1/ps < 
(AM S ) SM < 21.0/ps. 

We illustrate the correlations between AM S and BR(B° s d — > ^ + ^') in fig. 23 where we 
plot BR(B Q sd -> fj+fj,-) as functions of (AM s ) cxp /(AM s ) SM for tan/3 = 50 and M A = 300 
and 400 GeV by scanning the other MSSM parameters with the restriction that sparticles 
are heavier than 500 GeV and the B — > X s ^ constraint is satisfied. For each point in 
the MSSM parameter space we satisfy the experimental constraints on ex, V u b/V c b, a^x s 
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0.5 1 1.5 0.5 1 1.5 



AM S /AM S SM AM S /AM S SM 

Figure 23: Correlation between AM S /(AM S ) SM and B° d -> in the MSSM with 

flavour violation ruled by the CKM matrix. Lower (upper) branches of points correspond 
to < l + f s < 1 (1 + fs < 0). Current experimental bounds: BR{B Q S — ► < 2-10 -6 

(CDF) [35] and BR(B° d -> /i + /i") < 2.1 • KT 7 (BaBar) [34] are shown by the horizontal 
solid lines. Solid (dashed) vertical lines show the lower limit on AM S /(AM S ) SM following 
from eq. (7.62) with a = 0.52 as in ref. [23] (a = 0.71 as in [44]). 
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and AMd by appropriately adjusting the Wolfenstein and nonperturbative Bk, F Bd and £ 
parameters, that is, performing the unitarity triangle analysis [23, 7, 44] described in sec. 
6.4 (points for which these constraints cannot be satisfed are rejected). (AM d ) cxp and the 
parameter ek do not constrain our scan significantly as the corresponding quantities fd 
and f £ are small in most of the parameter space. Since the element V t e d s is well constrained 
by €r, Vub/Vcb, ct>4>Ksi fitting the experimental value of AM d practically fixes F Bd . For this 
reason, the correlation of BR(B d — > fi + with AM s /AMf p clearly seen in figures 23b 
and 23d is much tighter than the analogous correlation of BR(B° S — » (figures 23a, c) 

where the F Bs ~ £,F Bd can still be varied independently (some structure in the density 
of points visible in figs. 23a, c is an artifact of the scanning method used to produce the 
figure, preferring the points concentrated around the minimal or the maximal allowed 
value for BR(B Q S — > Vertical lines shown in figures 23 correspond to the lower 

limit on (AM s ) cxp /(AM s ) SM following from (7.62) with a taken from ref. [23] (solid) and 
[44] (dashed). 

For M A = 200 GeV and tan (3 = 50 shown in the plots of [29] all points corresponding 
to the rather unlikely scenario with 1 + f s < were eliminated by the combination of the 
lower limit (7.62) and the CDF upper bound BR(B° S -> //+//") < 2 x 10~ 6 [35]. This is 
not the case for heavier A and/or smaller tan (3 values. 

Therefore for such points we can only use (6.41) to find 



BR{B° d -> fi + fi') < 3.6 (3.1) • 10" 



1.15 



F Bs /F Bd 



BR(B° S -> /i+^ cxp 



10- 6 



(7.63) 



with the numerical factor corresponding to the analyses in [23] and [44], respectively. 
With the current CDF bound one has the upper bound BR(B d — > fi + fi~) < 8 (7) • 10~ 8 
which is still lower than the current BaBar bound [34]. 

For a more likely situation of < 1 + f s < 1 and (AM s ) exp satisfying (7.62) we get 
upper bounds on both branching ratios: 

BR(B° S ^ < 1.2 ■ 10~ 6 (8-10~ 7 ) for a = 0.52 (0.71), 
BR(B° d ^ < 3- 10~ 8 (2- 10" 8 ) for a = 0.52 (0.71). (7.64) 

where the two values for the parameter a correspond to the analyses in [23] and [44], 
respectively. This should be compared with the SM values that are in the ballpark of 
3 • 10 -9 and 1 • 10~ 10 , respectively. On the basis of our discussion of the contribution 
(AM S ) X± in sec. 6, we would like to emphasize that the upper limits on BR{B° sd — > fi + ^r) 
obtained here for heavy sparticle spectrum cannot be significantly altered by lowering the 
sparticle masses. 
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The upper bounds (7.64) and (7.63) are very impotrant because, as has been pointed 
out in ref. [7] and recently confirmed in [8], they do not apply in the MSSM in which 
flavour violation originates from the squark mass matrices. The reason is that in this case 
leading supersymmetric corrections to AM S and to BR(B® — > At + /i~) are governed by dif- 
ferent off diagonal elements of the squark mass matrices and the correlation between these 
two quantities does not hold. One concludes therefore, that within the supersymmetric 
framework, observation of the B® — > decay at the level higher than ss 3 x 1CT 8 , apart 

from implying that the scale of the Higgs boson sector is not far from the electroweak 
scale, would be also a very strong evidence for non-minimal flavour violation in the quark 
sector. 

In principle, if the measured values of BR(B® d — > turn out to be significantly 

bigger than predicted by the SM but do not respect the correlation with the measured 
value of AM S , this could also be a sign of nonminimal flavour violation. In practice, 
however, usefulness of this test is limited by the precision with which the nonperturbative 
parameters F Bs , B Bs etc. are known. 

8 Summary and Conclusions 

In this paper we have presented an effective Lagrangian formalism for the calculation 
of flavour changing neutral and charged currents that is suited for theories containing a 
sector of heavy fields. While our approach is rather general, we have discussed it explicitly 
in the context of the MSSM with the CKM matrix as the only source of flavour violation, 
heavy supersymmetric particles and neutral and charged Higgs masses below 500 GeV. 

Recently a number of analyses of large tan (3 effects in weak decays appeared in the 
literature. The present analysis is the first one that includes simultaneously 

• Resummed large tan (3 effects, 

• SU(2) x Z7 (1) gauge symmetry breaking effects, 

• The contributions of electroweak couplings g\ and g 2 to the effective parameters in 
the low energy Lagrangian. 

Our approach is therefore equivalent to the full diagramatic one loop calculation supple- 
mented by the resummation of the tan (3 enhanced terms. We have analyzed the impor- 
tance of these effects demonstrating that for reliable quantitative predictions all of them 
have to be taken into account. Thus, we have generalized previous analyses where either 
the resummed large tan/5 effects have been calculated in the SU{2) x U(l) symmetry 
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limit [20, 27, 28] or the complete one loop calculations have been performed leaving out, 
however, the resummation of large tan/5 effects [6, 22, 21]. 

To compare our approach with the one of [27, 28] we have also derived analytic formulae 
for the neutral and charged Higgs boson couplings to quarks in the SU (2) x U(l) symmetry 
limit resumming large tan (3 terms. Our results agree here with those presented in [27, 28]. 
These transparent formulae, which make all large tan/5 effects explicit, depend only on 
four universal flavour independent parameters e , ey, e' Q and e' Y . However, the comparison 
of the results obtained in this approximation with the ones following from the complete 
calculation shows clearly that the former can only be used for a semi-quantitative analysis. 
We have also shown that the accuracy of the approximation can be greatly improved 
while maintaining the transparency of the resulting expressions. This is achieved by 
introducing flavour dependence into the parameters eo, ey, e' and e' Y by extracting them 
directly from the complete (but relatively simple to compute) correction to the down 
quark mass matrix. In this manner we have obtained analytic formulae for the neutral 
and charged Higgs couplings to quarks that still make the large tan (3 effects explicit but 
include automatically the effects of electroweak couplings g 1: g 2 and of the, slightly less 
important, SU(2) x (7(1) symmetry breaking effects. Our analytic formulae reproduce 
within 5 — 10% the results of the full approach to be compared with only 20—40% accuracy 
of the formulae obtained in the approach of refs. [27, 28]. 

The main message of our analysis is that the simple approach adopted in [20, 27, 28] 
overestimates the effects of the resummation of large tan (3 contributions. This is clearly 
seen in figs. 18, 20. 

The most spectacular effects in weak decays predicted by the MSSM with large tan f3 
are the huge increase of the branching ratios BR(B® d — > and a significant sup- 

pression of AM S . As pointed out in [29] these effects are tightly correlated so that for 
most probable ranges of supersymmetric parameters the experimental lower bound on 
AM, implies an upper bound on BR(B Q sd — > fi + fi~). In [29] and in the present paper 
both these quantities, as well as their correlation have been calculated including for the 
first time all the effects listed above. Simultaneously we have clarified some points in the 
calculation of large tan f3 effects in the B — > X s 7 rate that we have used to constrain the 
allowed range of supersymmetric parameters. 

For the same set of parameters our results for BR(B® d — > without the inclusion 

of the AM S constraint are typically between the values presented in [27, 28] and in [6, 
21, 22]. We also find that the impact of the resummed large tan (3 corrections on our 
analysis of AM S in [23] is smaller than obtained in [27, 28]. Moreover we have shown that 
large tan (3 effects in AM, are dominated by double neutral Higgs penguin diagrams and 
receive significant contributions from box diagrams with charged Higgs exchanges. 
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For the near future of particular interest are the upper bounds on BR(^B^ d — > ji^ji ) 
as functions of the ratio (AM s ) exp /(AM s ) SM that we have already presented in [29]. In 
the most likely scenario with < (1 + f s ) < 1 these bounds are becoming very strong 
when this ratio approaches unity. This is illustrated in fig. 23. For (AM s ) cxp > (AM S ) SM 
substantial enhancements of BR(B Q sd — > with respect to the values obtained in 

the SM are not possible within the MSSM scenario considered here. Violation of our 
bounds would either rule out this version of MSSM and/or signal new sources of flavour 
violation [7]. 

As the upper bounds on BR(B° s d — > discussed here are sensitive functions of 

the ratio (AM s ) exp /(AM s ) SM , their quantitative usefulness will depend on the value of 
(AM s ) cxp and on the accuracy with which (AM S ) SM can be calculated. In this respect 
the present efforts of experimentalists to measure BR(B® d — > ^ + ^T) and AM S and of 
theorists to calculate Fs d3 and the parameters Bd yS are very important. 

Note added: Very recently a new analysis of resummed large tan f3 effects in FCNC transi- 
tions has been presented in [49]. It goes beyond the framework of minimal flavour and CP 
violation and discusses the ratio e' / e and CP asymmetries not considered in our paper as 
the new physics effects in these quantities in our scenario are tiny. However these authors, 
similarly to [27, 28], work in the SU(2) x U(l) limit and do not take into account the 
refinements discussed in section 5 nor they included the constraint from B — > X s 7 decay 
on the supersymmetric parameters. On the other hand they stress that the inclusion of 
light quark contributions in the formalism of [27, 28] could have an impact for certain 
very special values of supersymmetric parameters for which €3 approaches —1. We would 
like to emphasis that in our full approach of section 2 and in the improved approximation 
presented in section 5 the contributions of light quarks are automatically included. More- 
over when the experimental constraints from B — > X s 7 decay, AM Sj ^ and B — > /x 4 piT are 
taken into account and the Higgs masses are kept below 500 GeV, 63 never approaches 
— 1 and the light quarks except for assuring the GIM mechanism are unimportant in the 
resummed large tan (3 corrections. 
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Appendix Feynman rules and Green's functions 



For completeness we collect here some of the the Feynman rules given in [30] and formulae 
for the 2-, 3- and 4-point Green's functions used in the paper. 



A.l Feynman rules 

We write the fermion-fermion-scalar vertices as i(VlPl + VrPr)- The couplings Vl,Vr 
are listed below: 
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The three-scalar vertices are defined as iV, with V given by: 
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^ M W„ T/ J/ 7 (J+3)i* 7 (J+3)i* A lit 

-yujVdjV ^ D z u d 



t i 1 ( Z^l n, 7 7 lk I 7 lfc I 7- 7/ 2 7 2fc I v JI 7 Ij *7 Ji * 

+ - ^ - -D+ +(Z 1 H k ^y Uj V JI - Z^A^ J V KI )Z^Z^ 

+ (Z 1 H k A« I *V JK - Z%w dl V JI )Z J j*Z^* 

As explained in sec. 2, the CKM matrix V in the above rules is related to the "mea- 
sured" CKM matrix V e g as in (2.21) or, in the approximation of sec. 3, as in (3.25). The 
Yukawa couplings are related to the physical quark masses as: 

Vdi = 7z , - . 7j\7 Vuj = m uj (A.l) 

v d (1 + e/tan/3) t> u 

Definitions of the mixing matrices Zrj, Zu, Z± and Zjv can be found in [30]. Note that 
it is V, not V e ff that enters the squark mass squared matrices. It should be also remem- 
bered that computing Wilson coefficients of the higher dimension operators generated by 
charginos one should modify their vertices according to the rules formulated in sees. 2.6 
and 3.5. 

A. 2 Quark self-energies 

Definitions of the formfactors are given in fig. 1. Also, (S^) 7,7 = (£^ L ) J/ * 
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fc=l z=i 

A.3 Higgs-fermion effectve couplings 

The effective Higgs-fermion vertex is defined in fig. 3. The formulae given below are 
completely general and are valid also in the MSSM with flavour violation originating in 
the sfermion mass matrices. 

i) CP-even neutral Higgs boson vertex H^djdi. Using the general notation of 
eq. (2.24) with S = H% (#° = H°, H% = h°) we find: 

6 
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ii) CP-odd neutral Higgs boson vertex A^djdj. In this case S = A° k (A\ = A , 

A° 2 = G°) and 
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iii) Charged Higgs boson vertex H^ujdi where H + l = H + and = G + (see fig 4). 
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+ VSm Nl m Cm C {m 2 Un ,m 2 Nl , m 2 c J\ (A.9) 
A. 4 Box diagram contributions 

Box contributions to AF = 2 processes as given below are related to the Wilson coefficients 
Ci defined in eq. (6.4) through 

= G 2 F M^(V^V^) 2 Ci(M S usY) (A.10) 

i) Charged Higgs contributions (including the mixed iJ ± W^ =F diagram). In order to 
simplify the notation, we denote P[ R l = [Pj* R \ JI , P[r = [Plr\ ' U an d similarly for the RL 
couplings ([P H ], [P G ] are defined in eqs. (2.39), (2.40)). If not written down explicitly, 
the arguments of the loop integrals D and D 2 are (m 2 H+ ,m 2 H+ ,m 2 M ,ffi 2 N ). 
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Chargino contributions. Arguments of the loop integrals are (m^, , m^ fe , ) . 
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\ ^ \ ^ \rR.Jkn*\rR.Jlm*\rLIkm\rLIln rrri r-> 

-jL 2^ Kh/c Kh/c Kh/c v duc m c m rn Cn L>o 



i 6 2 

\ ^ \ ^ TrRJkn-kTrRJlrn-kTrLIkrriTrLIln^ r-> 

2^ ^dc/c ^h/c Kh/c v duc m c m m Cn -L>o 

k,l=l m,n=l 
i 6 2 

\ ^ \ ^ \rLJkn*\rLJlm*\/-RIkm\rRIln frri r) 

-7 L 2^ Kh/c Ka/c Kh/c Ka/c m Cm m Cn V 

k,l=l m,n=l 
I & 2 



16 



E\ 1 \rLJkn*\rLJlm*\rRIkm\ l rRIln^ r) 
2^ v duc v duc v duc v duc ™c m m Cn DQ 



k,l=l m,n=l 



\ ^ \ ^ \rLJkn-k\rRJlm*\7-LIkm\rRIln lrn r-> 

~o 2^ v duc v duc v duc v duc mc m m Cn Do 



6 2 



1 V""* V""* t / LJlm*T 7 RJkn-k-t r LIkm\ r Rlln r-> 

2 2^ 2^ K dl/C* ^dl/C* V dUC V dUC U 1 



(A.12) 



A. 5 Loop integrals 

We define three- and four-point loop integrals at vanishing external momenta as: 



(47T) 



1 r rl d h 

1 r* I 2 2 2\ / ex. /i. 



i k 2n 



(2n) d (k 2 - mi){k 2 - m 2 2 )(k 2 - m§) ' 

;D2 n {m 1 ,m 2 ,m 3 ,m i ) = 



(47r )2--V-l' ""2, - 3 , -4/ y (27r)d ( p _ m 2 )( p _ m 2 )(A; 2 _ m 2^ k 2 _ m 2) 

Explicit formulae for the C , C 2 , -D and _D 2 functions are as follows 



C (x,y,z) 



y 



log- + 



log - , 



(x-y)(z-y) x (x-z)(y-z) x 



(A.13) 
(A. 14) 

(A.15) 



x 



C 2 (x,y,z) = A + log — + 



log- + 



log - . (A. 16) 



/i 2 (x-y)(z-y) x (x-z)(y-z) x 

where A = + log 4^7^ — 1 and /i is the renormalization scale. In the flavour changing 
penguin diagrams A and /i dependence always cancels out after summation over squark, 
chargino and neutralino mixing matrices. 



D (x,y,z,t) 



y 



log 



y 



log 



(y - x)(y - z)(y - t) x (z - x)(z - y)(z - t) x 
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+ («-,)(« !,,)(« -.-) "*;■ (A - 17) 

D 2 (x,y,z,t) = -log- + - — log — 

(y - x){y - z){y - t) x {z - x){z - y){z - t) x 

t 2 t 

+ 7 77 77 7 log"- (A.18) 

(t-x)(t-y)(t-z) B x 
The two-point loop integrals B and B 1 at s — can be expressed as 

B o (0,x,y) = A + ^— log 4 + ^^ log A + log 4 + ^^ log- (A.19) 
£ — y fi z y — x /r [i £ x — y y 

B 1 (0, x, y) = l£b(0, y) - 1^ + ^-^^ log - = \ + \c 2 {x, y, y) (A.20) 
2 4x-|/ 2 (x — y) z y 4 2 

The function H 2 of eq. (5.4) is related to C as follows: 

H 2 (x, y) = -m 2 C (m 2 , a^ 2 , ym 2 ) = -C (l, x, y). (A.21) 
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